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RELATIVISTIC LINES AND REFLECTION FROM THE INNER ACCRETION DISKS AROUND NEUTRON STARS 
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ABSTRACT 

A number of neutron star low-mass X-ray binaries have recently been discovered to show broad, asymmetric 
Fe K emission lines in their X-ray spectra. These lines are generally thought to be the most prominent part 
of a reflection spectrum, originating in the inner part of the accretion disk where strong relativistic effects can 
broaden emission lines. We present a comprehensive, systematic analysis of Suzaku and XMM-Newton spectra 
of 10 neutron star low-mass X-ray binaries, all of which display broad Fe K emission lines. Of the 10 sources, 
4 are Z sources, 4 are atolls and 2 are accreting millisecond X-ray pulsars (also atolls). The Fe K lines are 
well fit by a relativistic line model for a Schwarzschild metric, and imply a narrow range of inner disk radii 
(6-15 GM/c 2 ) in most cases. This implies that the accretion disk extends close to the neutron star surface 
over a range of luminosities. Continuum modeling shows that for the majority of observations, a blackbody 
component (plausibly associated with the boundary layer) dominates the X-ray emission from 8-20 keV. 
Thus it appears likely that this spectral component produces the majority of the ionizing flux that illuminates 
the accretion disk. Therefore, we also fit the spectra with a blurred reflection model, wherein a blackbody 
component illuminates the disk. This model fits well in most cases, supporting the idea that the boundary layer 
is illuminating a geometrically thin disk. 

Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

X-ray emission lines from the innermost accretion disk 
are well-known in supermassive and stellar-mass black holes, 
where th ey are shape d by relativistic effects (for a re- 
view, see Millerl 120071) . The utility of such lines is two- 
fold: they can be used to constrain the spin of the black 
hole, and they can be used to constrain the nature of 
the innermost accretion flow, particularly the proximity of 
the disk to the blac k hole (B renneman & Reynolds 120061: 
Miniutti et al l 120091: iReis et al.1 l2009at iMiller et alJ 12 006a. 
2008b J2009tlRevnolds et alJ20 09: Wil kinson & Uttleyll2009l: 
Schmoll et al. 2009; Fabian et al. 20091). 



Disk lines are produced in a fairly simple manner: a source 
of hard X-ray emission that is external to the disk is all that 
is required. The specific nature of the hard X-ray emission 
- whether therm al or non-thermal, whe ther due to Compton- 
upscattering ( e.g.lGierlinski et alJI 19991) or emission from th e 
base of a jet dMarkoff & NowakH2004 iMarkoff et alj|2005l) . 
or even a hot blackbody - is less important than the simple 
fact of a continuum source with substantial ionizing flux. The 
most prominent line produced in this process is typically an 
Fe K line, due to its abundance and fluorescent yield, how- 
ever other lines can also be produced. The overall interaction 
is known as "disk reflection" and more subtle spectral fea- 
tures, including a "reflection hump " peaking between 20 - 
30 keV, are also expected (se e , e.g. iGe orge & Fabian M99U 
Magd ziarz & Zdziarskil 119951: iNavakshin & Kallmanl 120011: 
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iBallantvne et alJl200ltlRoss & Fabianll200l . 

Largely owing to a combination of improved in- 
strumentation and concerted observational efforts within 
the last two years, asymmetric Fe K disk lines have 
been observed in 8 neutron star low-mass X-ray bina- 
ries dBhattacharvva & StrohmaverJl2007l: ICackett et al 112008 : 
[ Pandel et al.l 120081: iD'Ai et alj I2009t ICackett et al.1 12 009a; 
PapittoetalJ 120091: iShaposhnikovet al.1 120091: IReis et ail 
2009b; lDiSalvo et al.ll2009h llaria et alj|2009h . In the case of 
neutron stars, disk lines can be used in much the same way 
as in black holes to determine the inner disk radius. The ra- 
dius of a neutron star is critical to understanding its equation 
of state; disk lines set an upper limit since the stellar surface 
(if not also a boundary layer) truncates the disk. In the case 
of neutron stars harboring pulsars, disk lines can be used to 
trace the radial extent of the disk a nd to obtain a magnetic 
field constraint ( Cackett et al.ll2009ah . Finally, in a number of 
neutron stars, hot quasi-blackbody emission (potentially from 
the boundary layer) provides most of the flux required to ion- 
ize iron; the fact that the disk intercepts this flux suggests that 
the inner disk is geometrically thin, or at least thinner than the 
vertical extent of the boundary layer. 

Disk line spectroscopy provides an independent view on the 
inner accretion flow in neutron star LMXBs, wherein some 
constraints have been derive d using X-ray timing (for a re- 
view, see Ivan der Klisl 120061) . The higher frequency part of 
so-called kHz QPO pairings ma y reflect the Keplerian orbital 
frequency in the inner disk (e.g. lStella & Vietri 1998), for in- 
stance. Trends in the coherence of the low frequency part of 
kHz QPO pairs may indicate the point at which the disk has 
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reached an inne r stable orbit (Barret et al. 2006), though see 
Mendez (2006) for an opposing view. Thus, X-ray timing 
and spectroscopy offer independent windows on the neutron 
star and inner accretion flow in LMXBs. It has been noted 
that with the assumption of Keplerian frequencies, combin- 
ing the upper kHz QPO frequency with the radiu s given by a 
line measurement yields a mass determination dPiraino et alj 
l2000t lCackettltlill2008h . 

In order to understand possible differences in the inner- 
most accretion flow onto neutron star low-mass X-ray bina- 
ries, and in order to derive joint constraints on fundamental 
neutron star parameters, it is important to analyze the full 
sample of disk lines using a consistent, systemati c approach. 
Herei n, we analyze archiv al data from Suzak u (Mitsud a et alj 
120071) and XMM-Newton dJansen et all 1200 lb . using consis- 
tent and systematic reduction and modeling procedures. Each 
spectrum is fit both with phenomenological models, and with 
physically-motivated disk reflection models. In addition to 
Fe K lines reported previously, we also present the discov- 
ery of broad lines in two additional sources, GX 17+2 and 
HETE J 1900. 1-2455. The total number of neutron star 
LMXBs with relativistic lines now stands at 10. 

2. SAMPLE OF NEUTRON STARS 

In this paper we analyze data from 10 neutron star 
LMXBs, eight of which h ave recently had detections o f broad 
Fe K emission fin e s {B hattachar vva & Strohmaverj 1 2007 1: 
Cackett et al l 120081: iPandel et all 120081: ID' Ai et all 120091 
Cackett etaLll2009at iPapitto et a"fll2009t IShaposhnikov et all 
2009; lReisetal.1 l2009bt iDi Salvo et alj 120091: llariaetafl 
2009). In addition, we also present evidence for relativistic 
Fe K emission lines in two further sources (GX 17+2 and 
HETE J1900.1-2455). The total sample consists of a va- 
riety of different types of neutron star LMXBs: four atolls 
(Ser X-l, 4U 1636-53, 4U 1705-44, 4U 1820-30), four Z 
sources (GX 17+2, GX 340+0, GX 349+2, Cyg X-2) and two 
accreting millisecond X-ray pulsars (SAX J1808.4— 3658, 
HETE J1900.1-2455, which are also atolls). The obser- 
vations we analyze here were obtained with either Suzaku 
or XMM-Newton, both of which have good effective area 
through the Fe K region (6.4 - 6.97 keV). In Table[T]we detail 
the observations presented here, as well as giving the distance 
to each source assumed throughout the paper. 

3. DATA REDUCTION 

Details of our data reduction for the Suzaku and XMM- 
Newton observations follows. We have attempted to follow a 
standard, common analysis procedure for all sources as much 
as possible. Where we have had to deviate from this, specific 
details for those observations are also given. 

3.1. Suzaku data reduction 

The data reduction was performed using HEASOFT v6.6.2, 
and the latest calibration files (as of June 2009). Suzaku con- 
sists of both a set of soft X-ray CCD detectors (XIS) as well as 
a separate hard X-ray detector (HXD). Details of the data re- 
duction for the XIS and HXD are given below. For all sources 
we adopted the following standard data reduction method, ex- 
cept where explicitly stated. Since ou r original analysis o f 
Ser X-l, GX 349+2 and 4U 1820-30 dCackett et al.ll2008l) . 
there has been significant advances in the analysis tools for 
Suzaku, for instance one can now generate response files for 
any source region (allowing for pile-up correction). 



3.1.1. XIS data reduction 

Firstly, we reprocessed the unfiltered event files using 
xispi which updates the PI values for the latest calibration. 
From there, we then created a cleaned, filtered event list using 
the standard screening criteria provided by Suzaku. In addi- 
tion to the standard screening we were also careful to filter 
out any periods where the telemetry was saturated. This only 
occurred in a very small number of cases, which are noted be- 
low. There are event lists for each data mode used on each 
detector. We treated all of them separately until after spectra 
and responses were generated. 

There are a couple of exceptions to this procedure. When 
a source is bright, the script to remove flickering pixels 
(cleansis) sometimes produces spurious results, rejecting 
all the brightest pixels in the image. This was the case for a 
portion of both the second Cyg X-2 observation and the SerX- 
1 observation. To overcome this, one must run cleansis 
without any iterations (as recommended in the cleansis user 
guide). We adopted this procedure for the few observations 
where this occurs. 

From the cleaned, filtered event files we extracted spec- 
tra. Source extraction regions were chosen to be a rectan- 
gular box of 250 x 400 px for 1/4 window mode observations 
and 125 x 400 px for 1/8 window mode observations. For 
the brightest sources where pile-up is present, we excluded 
a circular region centered on the source, the radius of which 
depended on severity of pile-up (radii are given below). 

After extraction of spectra, we generated rmf and arf 
files using xisrmfen and xissimarf gen (using 200,000 
simulated photons) for each spectrum. We combined all the 
spectra and responses for the separate data modes (e.g., 2x2, 
3x3) used for each detector using the addascaspec tool. 
Additionally, we then combined the spectra from all available 
front-illuminated detectors (XIS 0, 2 and 3). The spectra 
were then rebinned by a factor of 4 to more closely match the 
HWHM of the spectral resolution. 

We now give observation specific details for each source: 

Ser X-l: a central circular region of radius 30 px was ex- 
cluded from the extraction region. There is also a small 
amount of telemetry saturation. We exclude times when this 
occurs, which corresponds to removing 3% of XIS 0, 2 and 3 
events, and 0. 1 % of XIS 1 events. 

4U 1705—44: in the first and third Suzaku observations no 
pile-up correction was needed. However, in the second obser- 
vation, when the source was significantly brighter, a central 
circular region of radius 20 px was excluded. In the first ob- 
servation 2 type-I X-ray bursts were seen in the lightcurve, 
which we filtered out. No bursts were detected in the other 
two observations. 

4U 1820—30: a central circular region of radius 40 px 
was excluded from the extraction region. There is also some 
telemetry saturation, and we exclude times when this occurs. 
This removes approximately 16% of events on all XIS detec- 
tors. 

GX 17+2: a central circular region of radius 30 px was ex- 
cluded from the extraction region for both observations of this 
source. Only XIS and XIS 3 data are analyzed - XIS 2 was 
not in operation, and XIS 1 was not operated in window mode. 
We detected two X-ray bursts (one during each observation), 
which were excluded. 

GX 349+2: a central circular region of radius 30 px was 
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TABLE 1 
Observation details 



Source 




Class 


Mission 


Obs. ID 


Obs. start date 


Exp. time (ks) 


Mode 


Ref. 


Distance (kpc) 


Distance ref. 
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4U 1820-30 




A 


Suzaku 


401047010 


14/09/2006 


10/30 


1/4 W, 1.0s 


1 


7.6±0.4 


13 


GX 17+2 




Z 


Suzaku 


402050010 


19/09/2007 


5/15 


1/4 W, 0.5s 


5 


9.8 ±0.4 


12 








Suzaku 


402050020 


27/09/2007 


6/18 


1/4 W, 0.5s 


5 






GX 340+0 




z 


XMM-Newton 


0505950101 


02/09/2007 


33 


PN: timing 


6 


1 1 ± 3 


14 


GX 349+2 




z 


Suzaku 


400003010 


14/03/2006 


8/20 


1/8 W, 0.3s 


1 


5 


11 








Suzaku 
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19/03/2006 


8/24 


1/8 W, 0.3s 


1 












XMM-Newton 


0506110101 


19/03/2008 


10 


PN: timing 


7 






CygX-2 




z 


Suzaku 
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16/05/2006 


<0.5 


1/4 W, 1.0s 


8 


11±2 
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z 


Suzaku 
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1/4 W, 0.5s 


5 
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Suzaku 


903003010 


02/10/2008 


21/31 


1/4 W, 1.0s 


9 


3.5±0.1 
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30/09/2008 


43 
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9,10 






HETE J 1900.1 


-2455 
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Suzaku 


402016010 


16/10/2007 
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1/4 W, 2.0s 


5 


3.6±0.5 


12 
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NOTE. — A: atoll, Z: Z source, AMXP: accreting millisecond X-ray pulsar. The Suzaku exposure times given are (1) the good time for the individual XIS 
detectors, (2) the good time for the HXD/PIN. For Suzaku modes 'W 1 stands for window-mode (where a sub-array is used). The number following that is the 
frame time. The nominal frame time for 1/4 W is 2.0s, and for 1/8 W is 1.0s, shorter frame times encur deadtime. Where no uncertainty in distance is given, we 
assume a 25% uncertainty. 

excluded from the extraction region for both observations of 
this source. The XIS 0, 2 and 3 detectors were operated in 
1/8 window mode, with a 0.3s integration time. The XIS 1 
detector, however, was operated in 1/4 window mode, but also 
with a 0.3s integration time, and thus has a smaller livetime 
fraction (15% compared to 30%). 

Cyg X-2: the fir st observation of this source has previously 
been published by Shaposhnik ov et all (12009ft . These authors 
noted that a large fraction of this observation was affected 
by telemetry saturation and thus only used data from times 
when high telemetry rates were utilized. On re-analysis, we 
found that even during high telemetry rates the data still suf- 
fered telemetry saturation. Strictly filtering the data using 
time where there was no telemetry saturation leads to less than 
500s of good time for each detector. Due to the resulting low 
count statistics, we do not analyze this observation further. 

During the second observation, Cyg X-2 happened to be 
observed in a particularly bright phase, and we therefore ex- 
cluded a central circular regi on of r adius 75 px (see later dis- 
cussion of pile-up in Section I4.1.U . We only extract spectra 
from XIS 1 and XIS 3. XIS was operated in 'PSUM' (tim- 
ing) mode. This mode is not yet calibrated, thus the spectra 
cannot be used for analysis. There was significant telemetry 
saturation during this observation, and we remove all times 
when this occurs. This leads to the exclusion of 36% of events 
from XIS 1 and 9% from XIS 3. 

SAX J 1808.4-3658: we use the spectra from lCackett et alj 
(2009a), and do not reprocess the data here. 

HETE J 1900.1 —2455: we only use the sections of the data 
that were operated in 1/4 window mode. Given the low source 
count rate, no pile-up correction was required. 

3.1.2. HXD/PIN data reduction 

In every case we downloaded the latest 'tuned' background 
file from the Suzaku website. The PIN spectrum was extracted 



from the cleaned event file using good time intervals com- 
mon to both the PIN event file and the 'tuned' background 
file. The resulting source spectrum was then deadtime cor- 
rected. The non-X-ray background spectrum was extracted 
from the 'tuned' background file, using the same common 
good time intervals. In addition to the non-X-ray background, 
the cosmic X-ray background also contributes a small (~ 5%) 
amount to the total background. We fake a cosmic X-ray 
background spectrum following the standard Suzaku proce- 
dure, and add this to the non-X-ray background spectrum to 
create the total background spectrum. The response of the 
detector has changed over time, thus, we make sure that we 
used the correct response associated with the epoch of the ob- 
servation. In the few cases where type-I X-ray bursts were 
detected in the XIS lightcurves, we also filtered out the same 
times from the HXD/PIN data. 

3.2. XMM-Newton data reduction 

We processed the observation data files (ODFs) for each ob- 
servation using the XMM-Newton Science Analysis Software 
(v8). For the analysis here we only use timing mode data from 
the PN camera as it is the appropriate observing mode for the 
sources of interest. Calibrated event lists were created from 
the ODFs using the PN processing tool epproc. Exceptions 
to this procedure are the observations of SAX J1808.4— 3658 
and GX 340+0. We obtained the calibrated event lists for 
these observations directly from the XMM-Newton SOC. 

Before extracting spectra, we checked for periods of high 
background by extracting a lightcurve from an off-source strip 
with energy > 10 keV. We note where we found high back- 
ground levels below. Using only times of low background, 
we extracted spectra with the following criteria: quality flag 
= 0, pattern < 4 (singles and doubles) and in the energy range 
0.3 - 12 keV. The source extraction regions used were 24 pix- 
els wide (in RAWX), covered all RAWY values and were 
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centered on the brightest RAWX column. In a few cases, 
where the source was particularly bright, we had to exclude 
the brightest two RAWX columns (details noted below). The 
corresponding rmf and arf were created using the rmf gen and 
arf gen tools following the SAS procedures. 
Observation specific details follow: 

SerX-1: during the first observation, there is only a slightly 
high background at the beginning, which we choose to in- 
clude. There is no background flaring in the other two obser- 
vations. 

4U 1636—53: there were three X-ray bursts during the first 
observation, one during the second observation, and one dur- 
ing the third observation, all of which were excluded. Only 
during the second observation was there some background 
flaring. Removing this reduced the exposure time by 3.5 ks. 

4U 1705—44: one type-I X-ray burst occurred during 
this observation, and was excluded. We find no back- 
ground flaring . Note that this is not the dataset presented in 
iDi Salvo et alj (120091) . which is still proprietary at the time 
of writing. Other archival XMM-Newton observations of this 
source were observed in imaging mode, and so are not ana- 
lyzed here. 

GX 340+0: some low-level background flaring, which re- 
duces the exposure time by 7 ks. The source is particularly 
brigh t during this observation (see lightcurve in iD'Ai et alj 
2009), we therefore remove t he two brightest c olumns from 
the extract region. Given that D' Ai et aIJd2009h demonstrate 
that the Fe K line properties do not vary significantly with 
source state, we chose not to split the observation based on 
source state, but instead prefer to analyze the entire time av- 
eraged spectrum. 

GX 34 9+2: the source count rate is high during this ob- 
servation dlaria et al.ll2009t) . Following llaria et all d2009l) . we 
therefore exclude the two brightest columns from the extrac- 
tion region. There is no significant background flaring. An- 
other archival observation of this source is also available, 
however, it was observed in imaging mode, and thus is not 
analyzed here. 

SAX J1808.4-3658: we use the spectra from lCackett et alj 
(2009a) and do not reprocess the data here. 

4. SPECTRAL FITTING 

4.1. Phenomenological models 

Our first step is to fit the spectra with phenomen ological 
mode ls (throughout the paper we use XSPEC vl2; lArnaudl 
1996). There has been debate for many years as to the most 
app ropriate continuum model to fit to neu tron star LMXBs 
(e.g lWhite et al.lll988l:lMitsuda et alJl!989h . as it is often the 
case that sev eral different models fit the data equally we ll (e.g 
lBarrelll2001l) . A recent investigation bv lLin et al] d2007l) stud- 
ied two atoll sources over a wide range in luminosity. They 
investigated the luminosity dependence of measured tempera- 
tures for different continuum models, looking to see for which 
model choice the measured temperatures followed L T 4 . 
Based on this, they suggested a prescription for spectral fit- 
ting, using an absorbed disk blackbody (xspec model: diskbb 
lArnaudl["l996l) plus single temperature blackbody plus a (bro- 
ken) power-law for the soft and intermediate states, and a 
blackbody plus broken power-law for the hardest states. We 
have also foun d this model to fit the spe ctra of neutron star 
LMXBs well dCackett et al.l l2~008. 2009a!j). Here, we adopt 
this same model (though use a power-law, rather than a bro- 
ken power-law), testing that the addition of each component 



statistically improves the fit. 

We interpret the hot blackbody component as the boundary 
layer between the accretion disk and neutron star surfac e (e.g. 
Inoga mov & Sunyaevl[l99 9: Popham & Sunyaev 2001). One 
can also successfully model this component with a Comp- 
tonized model, such as comptt. When doing this, however, 
one gets high optical depths - the Comptonized component 
is close to a Wien spectrum. Rather than fit the more compli- 
cated Comptonization model, we chose the simpler blackbody 
model. 

The millisecond X-ray pulsar, HETE J 1900. 1-2455 has a 
much harder spectrum than the other sources we study here. 
In fact, we find that its broadband continuum spectrum can be 
well fit be a single power-law only. 

We found that several of the sources have an apparent emis- 
sion feature at around 1 keV, as has previously been observed 
in a number of neu tron star LMXBs, with a number of dif- 
ferent missions re.g. lVrtilekeTaT] [r988: Kuul kers et alJl!997h 
Schulz Hl999b . It has suggested that this f eature is a blend of a 
number of emission lines from Fe and O dVrtilek et alJ[l 988). 
Here, if the feature is present, we fit with a single Gaussian 
emission line. 

The presence of a broad Fe K emission line in eight of 
these sources has already been established elsewhere (see ref- 
erences earlier). Of the other two sources (GX 17+2 and 
HETE J1900. 1—2455) a Gaussian Fe K emission has pre- 
viously been reported in GX 17+2 dDi Salvo et all [2000b; 
iFarinelli et ail 12005b ICackett et afl l2009bh and the fact that 
a broad red-wing has not been previou sly observed is likely 
due to sensitivity (Cacke tt et al.ll200 9bj). Here we show that 
the profile is asymmetric. In HETE J1900. 1—2455, we de- 
tect a broad Fe K emission line for the first time, though the 
line can be equally well fit by a broad Gaussian or a relativis- 
tic line due to the low signal-to-noise ratio. Therefore, all 10 
objects studied here are seen to have broad Fe K lines. 

To fit these emission lines, we use t he relativistic line model 
for a Schwarzschild metric (diskline: lFabian et alj|1989l) . The 
Schwarzschild metric should be a good approximation for 
neutron stars - for neutron stars with reasonable equations of 
state and spin frequencies less than 600 Hz, the dimension- 
less angular momentum parameter is expected to be less than 
0.3, thus a ny deviations from the Schwarzschild metric are 
minor (see Mil ler et aLl[l9 98, for further discussion of this). 
We note that the diskline profile does not account for light 
bending, and so does not fully account for all relativistic ef- 
fects. However, it still remains the most appropriate available 
model here. 

We include this model in addition to the continuum model 
discussed above. The emission line parameters are the line 
energy, disk emissivity index, inner disk radius, outer disk ra- 
dius and inclination. We constrain the line energy to be within 
6.4 to 6.97 keV, the allowed range for Fe Koc emission. The 
disk emissivity index and inner disk radius are left as free pa- 
rameters in the fits. In most cases, the inclination is left as a 
free parameter (we discuss the exceptions to this later). Fi- 
nally, we fix the outer disk radius at 1000 GM/c 2 . Typical 
disk emissivities mean that the Fe K emission is quite cen- 
trally concentrated, and hence the outer disk radius is usually 
not well constrained, and has little affect on the measured in- 
ner disk radius. 

Several of the sources considered here have multiple obser- 
vations by the same, or different telescopes. The inclination 
of the source should not change over time, thus, we should ex- 
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pect to recover the same inclination from the Fe K line mod- 
eling for each observation. In order to investigate the extent 
to which this is true, we chose to let the inclination parameter 
vary from observation to observation. 

For the Suzaku spectra, we fit the XIS spectra over the en- 
ergy range 1.0 - 10 keV, ignoring 1.5 - 2.5 keV where large 
residuals are always seen (regardless of count rate), presum- 
ably of an instrumental nature. The PIN spectra are fit from 
12 keV upwards, with the upper energy depending on where 
the background dominates over the source spectrum. The 
combined front-illuminated spectrum, XIS 1 spectrum and 
HXD/PIN spectra are all fit simultaneously, with a floating 
constant used to account for any absolute flux calibration off- 
set between them. When fitting the XMM-Newton PN spectra, 
we fit over the 0.8 - 10 keV range. Source specific details, in- 
cluding exceptions to this standard procedure follow. 

SerX-1: SuzakufPIN spectrum fit from 12-25 keV. In the 
XMM-Newton spectra there are significant residuals around 1 
keV that are well fit by a Gaussian emission line. 

4U 1636—53: T he Fe line profile observed in this source 
(Pandel et al. 2008), is unlike what is seen in any of the other 
sources. It has a blue wing that extends to high energies, sug- 
gestive of a high inclination. Our spectral fitting does give a 
best fit wi th a high in c linatio n (we find that all cases peg at 
90°) and iPandel et alj J2008) also reported high inclinations 
(greater than 64° in all cases). However, no eclipses have 
been detected in this source and opt ical measurements are not 
especially constraining (36° - 74°. lCasares et alj|2006l) . 

4U 1705-44: Suzaku/PTN spectrum fit from 12 - 30 keV. 
The first Suzaku o bservation has qui te low sensitivity through 
the Fe K band (see lReis et alj200"9bT) . We therefore fix several 
of the line parameters. 

When fitting the XMM-Newton data of 4U 1705-44, we fit 
in the range 1.2 - 10 keV, as a relatively high Nh and low 
source flux lead to low statistics below 1.2 keV. We found 
that both the diskbb+bbody or bbody+power-law model fit the 
continuum spectra well, with the diskbb+bbody model having 
a marginally better % 2 value. The spectrum is of fairly low 
quality through the Fe K region, and thus a simple Gaussian 
models the Fe K line well. However, to allow comparison 
with other observations we still fit the line with the diskline 
model. When including the diskline model both choices of 
continuum model give consistent Fe K parameters, with al- 
most identical % 2 values: 1762.4 vs. 1762.5 (dof = 1752 in 
both cases) for diskbb+bbody and bbody+power-law contin- 
uum choices respectively. 

An unphysically high inclination is determined from the 
XMM-Newton observation. This source is not known to be 
eclipsing or dipping. There is a di fferent incl i nation deter- 
mine d in the Suzaku obs ervations dReis et alj l2009bl) than 
from iDi Salvo et al] (120091) . We revisit this issue later when 
discussing systematic uncertainties in Fe K modeling. 

4U 1820-30: SuzakufPIN spectrum fit from 12 - 25 
keV. With the improved calibration, the combined front- 
illuminated spectrum has a significan tly higher signal-to - 
noise ratio than the spectra presented by Cac kett et al.l d2008). 
With this improved spectrum, the detected emission line is 
found to be weaker than before, but still a significant detec- 
tion: with the emission line the best fit has % 2 = 1269 for 1 108 
degrees of freedom, whereas with no emission line the best fit 
has x 2 = 1362 for 1113 de grees of freedom. 

GX17+2: SuzakufPIN spectrum fit from 12 - 40 keV. 

GX 340+0: This source is hampered by a very large col- 



umn density (Nh ~ 10 23 cm -2 ). This means that even the 
neutral Fe K edge at ^7 keV is significant. Unfortunately, a 
large neutral Fe K edge can be a big problem when fitting the 
Fe K emission line. If there is any small change in ISM Fe 
abundance towards the source, it will change the depth of the 
edge, and consequently the appearance of the emission line. 
Additionally, the edge is likely to be more complicated than 
the simple step-like function used in absorption models. De- 
tailed structure has been observed for Fe L, and several Ne 
and O edge s with sensitive gratings spectra (e.g. lJuett et afl 
12004 120061) . Thus, with a large Nh, structure in the edge 
might affect results. 

When fitting these data, ID' Ai et all (|2009) chose to allow 
the ISM Fe abundance to vary, effectively changing the size 
of the Fe K edge. This produces a good fit, and the resulting 
Fe K emission line profile looks good. We, similarly, achieve 
good fits by varying the Fe abundance. Nevertheless, there is 
not a unique model that fits the data well (see the line pro- 
files in Fig.Q]). Fig.QJ shows the resulting line profile when 
the continuum is fit from 2.5 - 5 keV and 8-10 keV, ig- 
noring the Fe K region (and assuming standard abundances). 
Evidently, the resulting residuals do not look like a typical 
Fe K emission line p rofile, with a clear dip at around 7 keV. 
As mentioned above, iD'Ai et al] d2009t) successfully correct 
for this by assuming that the ISM Fe abundance towards this 
source is lower (see Fig. [TJ;). With a lower Fe abundance, 
the absorption edge is weaker, thus the line profile no longer 
shows a dip. 

Here, we alternatively suggest that this can also be in- 
terpreted as an Fe XXVI absorption line at 6.97 keV super- 
posed on a broad Fe K emission line (see Fig. QJ>, and Ta- 
ble^. There are p recedents for such a scenario. For instance, 
Boiri n et alj d2005l) find Fe XXV and Fe XXVI absorption lines 
superposed on a broad Fe K emission line in the dipping 
LMXB 4U 1323—62, and recently similar absorption lines su- 
perposed on a broad Fe K emission line were observed in the 
black hole binary GRO J1655-40 dDiaz Trigo et al.ll2007l) . 

As is apparent in Fig. [T] the two models lead to differ- 
ing emission line profiles. While in the model varying Fe 
abundance the line profile looks asymmetric, it is not as ob- 
viously asymmetric when including an absorption line. Fit- 
ting the emission line profiles with the diskline model we 
also achieve different measure ments of R^: 9.3 ± 0.2 GM/c 2 
(consistent with the analysis of lD'Ai et al.ll2009l) compared to 
22 ±4 GM/c 2 for the alternative model with an absorption 
line included. Due to the ambiguous nature of the line profile 
in this source we do not consider this spectrum further here, 
though parameters from the above spectral fits can be found 
in the Table [2] 

GX 349+2: We fit the two Suzaku spectra separately, even 
though the observation times are close together, in order to 
investigate any variability. The SuzakufPIN spectrum was fit 
from 12-25 keV. In the XMM-Newton spectrum, there are 
large resi duals at ~ 1 keV, that are well fit by a Gaussian emis- 
sion line dlaria et alj|2009l) . 

Cyg X-2: XIS spectra fitted from 0.8 - 10 keV, ignor- 
ing 1.5 - 2.5 keV. The SuzakufPIN spectrum was fit from 
12-25 keV. We find large residuals at ~ 1 keV (in all 
XIS spectra), that are well fit by a single Gaussian emis- 
sion line. This spect ral feature has prev i ously been observed 
in this source (e.g IVrtilek et alj 119881: IPiraino et al] 120021: 
IShaposhnikov et alj|2009l). and may be a complex of blended 
emission lines dVrtilek et aT] [T988: Sch ulz et al.l l2009). 
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TABLE 2 

Spectral fitting parameters for GX 340+0 



Parameter 


Variable Fe abundance 


Fe XXVI Absorption line 


N H (10 22 cm" 2 ) 


10.9 ±0.1 


10.0 ±0.1 


Fe abundance (rel. to solar) 


u - z:) -0.07 


1.0 (fixed) 


Disk kT (keV) 


1.64 ±0.04 


1.21 ±0.02 


Disk normalization 


127 ±7 


325 ±9 




2 44 + n 00 


1 93 ±0 01 


Bbody normalization (10~ 2 ) 


9.2±0.5 


15.2±().l 


Gaussian line energy (keV) 




6.97 (frozen) 


Gaussian a (keV) 




0.05 (frozen) 


Gaussian normalization 




(-1.1 ±0.1) x 10" 3 


Diskline line energy (keV) 


6.97_ .oi 


6.89±0.01 


Emissivity, q 
R m {GM/c 2 ) 


3.7±0.1 


1.85±0.14 


9.3±0.2 


22±4 


if) 


22±1 


37±4 


Diskline normalization (10 3 ) 


13.1 ±0.2 


5.5 ±0.2 


EW (eV) 


80 


38 


X 2 /dof 


1827/1548 


1925/1548 



NOTE. — The outer disk radius was fixed at 1000 GM/c~ in both cases. Normaliza- 
tion of the disk blackbody is in fi 2 nner cos//(Z)10) 2 where R\ nncT is the inner disk radius 
in km, D10 is the distance to the source in units of 10 kpc, and i is the disk inclination. 
Normalization of the blackbody component is L39/(D10) 2 , where L39 is the luminosity 
in units of 10 39 erg s -1 . The normalization of both the Gaussian and the diskline is the 
line flux in photons cm -2 , s , 

SAX Jl 808.4-3658: we do not re-fit the Suzaku and XMM- 
Newton spectra of th is source here, and refer the reader to 
iCackettet al.l d2009a) regarding specifics. 

HETE J 1900.1 -2455: The Suzaku/PIN spectrum was fit 12 
- 30 keV. This spectrum is quite noisy, though an Fe K emis- 
sion line is apparent. A broad (a = 0.8 keV) Gaussian fits the 
line well (% 2 = 1 148, d.o.f. = 1228), and a relativistic disk line 
is not statistically required. We do, however, fit the diskline 
to allow comparison with other data, and it gives a similarly 
good fit (x 2 = 1 146, d.o.f. = 1226). Also note that the contin- 
uum model for this source only requires a power-law compo- 
nent and no disk blackbody, or blackbody are required. 

4.1.1. Investigating the effects of pile-up 

Of all the sources observed, Cyg X-2 had the highest count 
rate observed by Suzaku, leading to severe pile-up. The con- 
tinuum sha pe changes d ue to pile-up, becoming artificially 
harder (e.g. Davis 2001). In order to investigate the robust- 
ness of the Fe K line profile in Cyg X-2, we extracted spectra 
from multiple different regions. In all cases, we used a box 
250 x 400 px, but had a circular exclusion region of varying 
radius centered on the source. The radius of this circle ranged 
from px to 90 px. We fitted the spectra from 0.8 - 10 keV, 
ignoring the 1.5-2.5 ke V energy range, as usual. The spectra 
were allowed to have completely different continuum models, 
fit using an absorbed disk blackbody plus a single temperature 
blackbody. Additionally, we had to include a narrow Gaussian 
at 1 keV. Here, we do not fit the PIN data, thus no power-law 
component is required. We do not fit the PIN data as pile-up 
significantly effects the shape of the XIS spectrum, hence the 
broadband spectrum will not fit well unless the XIS spectra 
are corrected for pile-up. 

We found that even though each spectrum had different con- 
tinuum shapes (due to pile-up effects), the Fe K line profile 
was remarkably robust. This is clearly demonstrated in Tab. [3] 
and Fig. 12 In Tab.f3]we give the continuum and Fe K line pa- 
rameters from fitting the spectra with different size exclusion 
regions. Fig.f2]shows the resulting line profiles for each spec- 
trum. The Fe K line is both quantitatively and qualitatively 
consistent, regardless of the pile-up correction. Thus, while 
with pile-up the continuum may vary significantly from the 




Energy (keV) 

FIG. 1. — Ratio of data to continuum model for GX 340+0. Each panel 
shows the same spectrum, fit with a different model (see Tablef2]for spectral 
fits), (a) The continuum is fit from 2.5-5 keV and 8-10 keV, ignoring the Fe K 
region, with an absorbed blackbody + disk blackbody model. A clear double- 
peaked structure is apparent, (b) An absorption line at 6.97 keV is added to 
the model along with a relativistic Fe K emission line. The ratio shown is 
that of the data to the continuum + absorption line model, thus, showing the 
emission line profile only. The red, sol id line shows the emission line model, 
(c) Here we follow D'Aietal. (2009), allowing the ISM Fe abundance to 
vary, changing the depth of the Fe K edge from ISM absorption. The emission 
line model (red, solid line) is shown. 

'true' (not piled-up) continuum shape (in the most extreme 
cases it is obvious that the continuum model parameters are 
not realistic), it appears that the Fe K line shape can be ro- 
bustly recovered as long as the continuum can be fit well. To 
further demonstrate this, we also show the Fe K profile from 
the Suzaku observation of Ser X-l (Fig. |2). We again see that 
the line profile remains consistent as the radius of the circular 
exclusion region is varied from px to 90 px. 

4.1.2. Results from phenomenological fitting 

The spectral parameters from this phenomenological mod- 
eling are given in Tables |4] and [5] (the tables have been split 
into continuum and line parameters). Table|6]gives the 0.5 - 
25 keV unabsorbed fluxes and luminosities. A summary of all 
the Fe K lines observed in neutron star LMXBs can be seen 
in Fig. f3] 

The distribution of measured inner disk radii is shown in 
Fig.|4] We only use one measurement per source, choosing the 
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TABLE 3 

Investigating pile-up in Cyg X-2: spectral parameters for different extraction 

REGIONS 



Parameter 




Excluded circle radius 






px 


30 px 


60 px 


90 px 


N H (10 22 cm" 2 ) 


0.1 47 ±0.003 


0.136 ±0.004 


0.128 ±0.004 


0.130±0.005 


Disk kT (keV) 


1.595 ±0.006 


1.555 ±0.002 


1.52 ±0.02 


1.45 ±0.04 


Disk normalization 


124 ±2 


151 ± 1 


165 ±3 


190 ±20 


Blackbody kT (keV) 


9±1 


4.1 ±0.1 


2.9±0.2 


2.2±0.2 


Bbody normalization (10 ) 


41±2 


7.0 ±0.1 


4.2±0.1 


4.4±0.6 


Line Energy (keV) 


6.97-o.o2 


6.97-0.02 


6.97-o.o2 


6.97_o.o2 


Emissivity, q 


5.8 ±0.6 


5.7±0.5 


6.3 ±0.8 


6.5±\<> 


R in (GM/c 1 ) 


8.0±0.2 


7.8 ±0.2 


7.9±0.3 


8.0 ±0.5 


iO 


24.0 ±0.5 


24.3 ±0.6 


23.9 ±0.6 


23±1 


Line normalization (10~ 3 ) 


7.2 ±0.4 


8.4±0.4 


8.7±0.6 


8.4 ±1.0 


EW (eV) 


69 


76 


79 


74 


NOTE. — Spectral parameters 


from fitting the 


XIS 3 data. While the continuum 


is variable due 



to pile-up effects, the Fe K line parameters are remain consistent. Normalization of the disk black- 
body is in « 2 nnCT cos;/(D10) 2 where R inner is the inner disk radius in km, D10 is the distance to the 
source in units of 10 kpc, and i is the disk inclination. Normalization of the blackbody component is 
L39/(D10) 2 , where L39 is the luminosity in units of 10 39 erg s _1 . The normalization of the diskline 
is the line flux in photons cm -2 . s~ ! . 



1.05 




Energy (keV) 



1.1 



1.05 




Energy (keV) 



FIG. 2. — Left: Fe line profile for Cyg X-2, excluding a central circular region of varying radius. Black: exclusion circle radius = px, red: radius = 30 px, 
green: radius = 60 px, blue: radius = 90 px. The Fe K line profile remains consistent regardless of the extraction region used. Right Fe line profile for Ser X-l 
from the Suzaku observation. Only data from XIS 3 is shown. Colors denote the same size extraction regions as for Cyg X-2. 




4U 1636-53 




Energy (keV) 

FIG. 3. — A summary of the Fe K emission lines in neutron star LMXBs. Plotted is the ratio of the data to the continuum model. The solid line shows the 
best-fitting diskline model. Data from Suzaku is shown in black (combined front-illuminated detectors), data from XMM-Newton is shown in red (PN camera). 
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TABLE 4 

Phenomenological model parameters: continuum 



Source Mission Obs. ID N H Disk Blackbody Blackbody Power-law Gaussian 









10 22 cm" 2 


kT in (keV) 


Norm. 


kT (keV) 


Nomr. (10~ 2 ) 


r 


Norm. 


Line energy (keVl 


a (keV) 


Norm. (10 -2 ) 


Serpens X 1 


Suwku 


40 M 14 MHO 


0.67 ± 0.04 


1 .26 ± 0.01 


97 ±5 


2.27 ±0.02 


5.23 ±0.06 


3.4±0.3 


0.9 ±0.2 










XMM-Newton 


0084020401 


0.35 ±0.01 


0.95 ± 0.01 




1 .76 ±0.01 


4.46 ±0.01 






1 . 10 ±0.01 


0. 1 1 ± 0.01 








(Ml841P0^O 1 


35 ±0 01 


99 ±0 01 


17^ ~ ' 


1 8 1 ±0 01 


3 13 ±0 02 






1 10 ±0 01 


14 ± 01 


1 4 too' 




XMM -Newton 


0084020601 


0.35±0.01 


0.95 ±0.01 


203 ±2 


1.76±0.01 


3.80 ±0.01 






1.10±0.01 


0.13±0.01 


1 .2 ± 0.01 


4U 1636-53 


XMM-Newton 


0303250201 


0.36±0.01 


0.21 ±0.01 


(1.0±0.2) x 10 4 


1.97 ±0.03 


0.15 ±0.01 


1.9±0.1 


0.22±0.01 










XMM-Newton 


0500350301 


0.51 ±0.01 


0.89 ±0.01 


101 ±4 


2.00±0.01 


1.19±0.01 


3.8±0.1 


0.54 ±0.01 










XMM-Newton 


0500350401 


0.50±0.01 


0.96 ±0.01 


103 ±4 


2.03 ±0.01 


1.57±0.01 


3.8±0.1 


0.58 ±0.01 








4U 1705-44 


Snzukii 


401046010 


2.0 


0.09 ±0.01 


(2.3+[ 2 )xl0 7 


1.05 ±0.05 


0.063 ±0.005 


1.65 ±0.01 


0.124±0.002 










Sazxika 


401046020 


1.95 ±0.05 


1.16±0.02 


56±2 


2.12±0.03 


2.40 ±0.02 


2.85 ±0.08 


0.64 ±0.06 










Sazxika 


401046030 


2.02 ±0.05 


0.82 ±0.02 


82 ±5 


1.95 ±0.03 


0.95 ±0.01 


3. 00 ±0.06 


0.49 ±0.03 










XMM-Newton 


0402300201 


1.53 ±0.01 


1.44 ±0.03 


2.5±0.2 


3.2 ±0.3 


0.25 ±0.01 












4U 1820-30 




401047010 


0.21 ±0.01 


1.18±0.01 


103 ±5 


2.44±0.01 


7.7±0.1 


2.3±0.1 


0.37 ±0.03 








GX 17+2 


Sazxika 


402050010 


2.2 ±0.1 


1.75 ±0.02 


86±5 


2.66 ±0.02 


6.1 ±0.2 


2.6 ±0.6 


0.4 ±0.4 










Sitzxikit 


402050020 


2.14±0.05 


1.74 ±0.02 


92 ±4 


2.55 ±0.02 


9.2 ±0.3 


2.4 ±0.3 


0.2 ±0.2 








GX 349+2 


Sitzxikit 


400003010 


0.77 ±0.01 


1.53 ±0.02 


112±4 


2.24 ±0.02 


12.5±0.2 


1.3 ±0.7 


00 ,+0.05 
UUI -0.01 












400003020 


0.85 ±0.02 


1.48 ±0.02 


105 ±4 


2.30±0.OI 


9.6±0.1 


2.3±0.1 


0.31 ±0.01 










XMM-Newton 


0506110101 


0.67 ±0.01 


1.17 ±0.01 


258±6 


1.85 ±0.01 


16.7 ±0.3 






1.08±0.01 


0.09 ±0.01 


3.9 ±0.4 


Cyg X-2 




403063010 


0.20 ±0.02 


1.50±0.01 


163 ±4 


2.21 ±0.02 


3.6±0.1 


2.8±0.1 


60 ±20 


1.01 ±0.01 


0.11 ±0.01 


9.0 ±0.6 


SAX J1808.4-3658 


Suzuki! 


903003010 


0.046±0.04 


0.48 ±0.01 


5.9±0.3 


1.00 ±0.02 


0.127±0.001 


1.93±0.01 


0.20±0.01 










XMM-Newton 


0560180601 


0.23 ±0.01 


0.23 ±0.01 


?0400+ 13200 


0.40±0.01 


0.125 ±0.004 


2.08±0.01 


0.33 ±0.01 








HETE J 1900. 1-2455 




402016010 


0.032 ±0.04 








2.14±0.01 


0.063 ±0.001 









Note. — Continuum spectral filling parameters from phenomenologieal modeling. All uncertainties are lo. Where no uncertainly is given, thai parameter was fixed. The X~ values for the fits are those in Tab.[5] Normalization or the disk blackbody 
is in/? 2 iiner cosi/(Z)10) 2 where ff mner is the inner disk radius in km, D10 is the distance to the source in units of 10 kpc, and i is the disk inclination. Normalization of the blackbody component is Z.39/(Z)10) 2 , where L39 is the luminosity in units of 
10 39 erg s — 1 . Normalization of the power-law component gives the photons keV - 1 cm -2 s"" 1 at 1 keV. Finally, the normalization of the Gaussian is the total number of photons cm -2 s — 1 in the line. 



TABLE 5 

Phenomenological model parameters: Fe K line 



Source 


Mission 


Obs. ID 


Line Energy (keV) 


Emissivity index 


Diskline 
R m (GM/c 2 ) 


Inclination (°) 


Nomr. (10 -3 ) 


EW (eV) 


X 2 /dof 


Serpens X- 1 




401048010 


f>-97_ .02 


4.8 ±0.3 


8.0 ±0.3 


24 ±1 


5.6 ±0.4 


98 


1442/1105 




XMM-Newton 


0084020401 


6.66 ±0.02 


1.8±0.2 


25 ±8 


32 ±3 


1 .8 ± 0. 1 


43 


2125/1831 




XMM-Newton 


0084020501 




3.7±4 


14±1 


13 ±2 


1.6 ±0.1 


50 


1961/1831 




XMM-Newton 


0084020601 


6.66 ±0.02 


2.0±0.2 


26±8 


27 ±2 


1.4 ±0.1 


38 


2151/1831 


4U 1636-53 


XMM-Newton 


0303250201 


6.4+ - 02 


2.7±0.1 


6.1+0-4 


> 80 


1 .2 ±0.9 


184 


1939/1832 




XMM-Newton 


0500350301 


6 4+0.04 


2.8 ±0.1 


7.7 ±0.6 


>77 


1.2 ±0.1 


106 


2112/1832 




XMM-Newton 


0500350401 


6.4+0.02 


2.8±01 


6.0+ - 2 


> 84 


2.5±0.1 


164 


2406/1832 


4U 1705-44 


Suzaku 


401046010 


6.62 + 0.03 


3.0 


10.0 


24 


<0.44 


<76 


953/1122 




Suzahi 


401046020 


6.94 + 0.06 


4.8±0.6 


7.3 ±0.4 


24±4 


3.3±0.3 


120 


1266/1120 






401046030 


6.9 + 0.1 


3.7 ±0.3 


6+0.2 


15+3 


1.3 ±0,1 


94 


1152/1120 




XMM-Newton 


0402300201 


6.43 ± 0.03 


2.1±0.1 


6+ 1 


>74 


0.34±0.05 


151 


1760/1752 


4U 1820-30 




401047010 


6-97-0.03 


3.3 ±0.1 


6.0+°- 1 


< 8 


2.6±0.3 


36 


1269/1108 


GX 17+2 


Suznkti 


402050010 


6.8±0.1 


2.6±0.2 


7±3 


15±15 


6.0 ±0.5 


46 


636/606 




Suzxikii 


402050020 


6.75 + 0.07 


4.6±0.7 


8.0 ±0.4 


27 ±2 


14±1 


81 


628/606 


GX 349+2 


Suzaku 


400003010 


6-97_ .01 


3.9±0.2 


7.5±0.4 


26 ±1 


16 ± 1 


105 


1268/1106 




Suzxikti 


400003020 


6.64 + 0.03 


2.3 ±0. 1 


10±2 


42 ±2 


9.2±0.6 


79 


1262/1106 




XMM-Newton 


0506110101 


6.76 + 0.01 


1.9±0.1 


6.0+0- 5 


55±2 


12.1 ±0.3 


79 


2531/1831 


Cyg X-2 




403063010 


6-97_ .02 


6.5 ± 1.3 


8.1 ±0.9 


25 ±3 


7.0 ±0.6 


65 


1305/1129 


SAX J1808.4— 3658 


Suzaku 


903003010 


6.40+0.04 


2.9±02 


n±\ 


50^ 


0.70 ±0.06 


134 


2987/2550 




XMM-Newton 


0560180601 


6.4+0.06 


3.0±O2 


13±4 


59±4 


0.68 ±0.04 


118 


2124/1950 


HETE J 1900. 1-2455 




402016010 


6 - 97 -0.2 


5 + 1 


6+ 1 


20±4 


016 ±0.02 


122 


1146/1226 



Note. — Fe K line parameters from phenomenological modeling. All uncertainties are la. Where no uncertainty is given, that parameter was fixed. In the diskline model, the emissivity of the 
disk follows a power-law radial dependence, R~ a , where q is the emissivity index we report here. The normalization of the diskline is the line flux in photons cm -2 s — 1 . 



TABLE 6 

Source fluxes and luminosities in the 0.5 - 25 keV range 



Source 


Mission 


Obs. ID 


Phenomenological 


Reflectio 










Unabs. Flux 


Luminosity 


Unabs. Flux 


Luminosity 








(10 -8 ergs" 1 cm -2 ) 


(10 37 ergs" 1 ) 


(10 -8 ergs -1 cm -2 ) 


(10 37 ergs" 1 ) 


Serpens X-l 




401048010 


1.19±0.01 


10,0 + 5.0 


1.32 + 0.08 


11.1 ±5.6 




XMM-Newton 


0084020401 


0.70±0.01 


5.9 + 3.0 


0.71+0.01 


6.0±3.0 




XMM-Newton 


0084020501 


0.59±0.01 


5.0 + 2.5 


0.60 + 0.01 


5.1 ±2.5 




XMM-Newton 


0084020601 


0.65 ±0.01 


5.5 + 2.7 


0.66 + 0.01 


5.6±2.8 


4U 1636-53 


XMM-Newton 


0303250201 


0.17+0-01 


0.73 + 0.05 


0.18 + 0.02 


0.78 ±0.09 




XMM-Newton 


0500350301 


0.39±0.01 


1.7 + 0.1 


0.37 + 0.01 


1.6+0,1 




XMM-Newton 


0500350401 


0.48±0.01 


2.1 ±0.1 


0.46 + 0.03 


2.0 + 0.1 


4U 1705-44 


Suzaku 


401046010 


0.17+0-01 


0.68 ±0.06 


0.57+0.05 


2.3 + 0.3 




Suzaku 


401046020 


0.61 ±0.01 


2.5±0.2 


0.49 + 0.01 


2.0 + 0.1 




Suzaku 


401046030 


0.30±0.01 


1.2 ±0.1 


0.31+0.01 


1.2+0-1 




XMM-Newton 


0402300201 


0.043 ±0.001 


0.17 ±0.01 


0.042 ±0.001 


0.17 + 0.01 


4U 1820-30 


Suzaku 


401047010 


1.21 ±0.01 


8.4±0.9 


1.21 ±0.13 


8.4+1.3 


GX 17+2 


Suzaku 


402050010 


2.37±0.01 


27.2 ±2.2 


2.27 ±0.01 


26.1+2.1 




Suzaku 


402050020 


2.60±0.01 


29.9 ±2.4 


2.62 ±0.01 


30.1+2.5 


GX 349+2 


Suzaku 


400003010 


2.35 ±0.03 


7.0±3.5 


2.36 ±0.13 


7.1 ±3.6 




Suzaku 


400003020 


1.97 ±0.01 


5.9±2.9 


1.92+0.01 


5.7±2.9 




XMM-Newton 


0506110101 


2.38±0.01 


7.1 ±3.6 


2.40 ±0.09 


7.2±3.6 


Cyg X-2 


Suzaku 


403063010 


2.22 ±0.01 


32.1 ±11.7 


2.30±0.01 


33,3 + 12.1 


SAX J1808.4-3658 


Suzaku 


903003010 


0.20±0.01 


0.29 ±0.02 


0.27 ±0.02 


0.40 + 0.04 




XMM-Newton 


0560180601 


0.25 ±0.01 


0.37 ±0.03 


0.26 ±0.01 


0.38 + 0.03 


HETE J1900. 1-2455 




402016010 


0.034±0.001 


0.053 + 0.015 


0.032±0.005 


0.050 ±0.016 



Note. — Distances used arc given in TablefT] 
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FIG. 4. — Measured inner disk radii (GM/c 2 ) from phenomenological fits 
(left) and reflection fits (right). Where there are multiple observations, the 
one with the smallest fractional uncertainty was chosen. The upper x-axis 
shows the inner disk radius in km assuming a 1 .4 M© neutron star. The hori- 
zontal dotted lines separate the Z sources (top), atolls (middle) and accreting 
millisecond X-ray pulsars (bottom). The vertical dashed line marks 6 GM/c 2 , 
the lower limit on the inner disk radius allowed by the model. 

one with the smallest fractional uncertainty. The vast majority 
of the measured inner disk radii fall into a relatively narrow 
range of 6 - 15 GM/c 2 . Only two observations (of Ser X-1) 
fall outside this range, and even then the inner disk radii from 
those observations is only a little more than la away. Note 
that there are a number of fits in which R[ n pegs at the lower 
limit of the model (6 GM/c 2 , the innermost stable circular 
orbit in the Schwarzschild metric). This corresponds to about 
12 km, assuming a 1.4 M e neutron star. Given that this lower 
limit is very similar to predicted neutron star radii, smaller 
inner disk radii are not expected. The lines in the observations 
with R{ n = 6 GM/c 2 are well fit by the model, and do not 
show large residuals in the red wing of the line (which would 
be expected if Ri n were smaller than measured). 

One possibility for the large number of sources that reach 
the lower limit of the model is that there is some contribu- 
tion to the broade ning of the line from Compton scattering. 
iReis et all d2009bl) found that for 4U 1705-44 the measured 
inner disk radius increased when such effects were taken into 
account using reflection modeling. 

We note that the differences in equivalent width (EW) be- 
tween our original Suzaku analysis of Ser X-1, GX 349+2 and 
4U 1820—30 and this analysis can be attributed to updated 
responses. During the original analysis the standard response 
matrices had to be used, which did not fully account for the 
window mode, and therefore gave incorrect fluxes. 

4.2. Reflection models 

As noted earlier, the Fe Koc emission line is only the most 
prominent feature expected from disk reflection where hard 
ionizing flux irradiates the accretion disk leading to line emis- 
sion. In order to more self-consistently model the spectra, we 
fit reflection models which include line emission from impor- 
tant astrophysical metals (the strongest line being Fe Koc), as 
well as effects such as Compton broadening. 

As is clear from the phenomenological mo deling of the 
spect r a (and as has b een noted before, e.g., Cack eTtet al.l 
I200ailariaet!r]l2009l) . in the atoll and Z sources, the black- 
body component (potentially associated with the boundary 
layer) is often dominating between 8-20 keV. Thus, it is 
this emission that is providing the vast majority of ioniz- 
ing flux that illuminates the accretion disk forming the Fe 
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FIG. 5. — Blurred reflection models for illumination by a power-law, 
spectral index = 2 (black) and for illumination by a blackbody, kT = 2 
keV (red) both for an ionization param eter of 1000 (models used are from 
Ballantyne et al. 2001; Ballantyne 2004). The models are normalized so that 
they both peak at 1 . The spectra are blurred using the rdblur model (a con- 
volution with the diskline model), both assuming an emissivity index q = 3, 
inner disk radius = 10 GM/c 2 , outer disk radius = 1000 GM/c 2 , and inclina- 
tion = 30° . Both spectra are also absorbed assuming a quite typical value of 
N H = 5 x 10 21 cm" 2 . 

K line (and reflection spectrum). In order to model the re- 
flection component, one therefore needs to employ a reflec- 
tion model whe re a blackbody component provides the il- 
luminating flux (Ballantvne & Strohmaver 2004; Ballantyne 
2004). This differs from the stand ard reflection models 
used to model black hole sources (e.g. Ball antvne et al J|200 It 
Ross & Fabian 2007) which assume illumination by a power- 
law spectrum. The r esulting reflection spectra are quite dif- 
ferent (see FigfS] and Ballantvne 2004, for details), with the 
blackbody reflection model dropping off quickly above 10 
keV compared to the power-law reflection. 

Here, we use models created to study the Fe K emis- 
sion line observed in 4U 1820—30 during a superburst, 
where a constant-density slab is i lluminated by a b lackbod>0 
dBallantvne & Strohmavel 12004 iBaflantvnd [20041) . We use 
the model assuming solar abundances for all sources except 
4U 1820—30, where we use t he model for a He system 
dBallantvne & Strohmaverll2004 . We note here that the den- 
sity used in these models is «h = 10 15 cm~ 3 . This is signif- 
icantly lower than th e density expected at the inner disk in 
neutron star systems. Ballantvne (2004) tested increasing the 
density and found that it did not significantly change the Fe 
K emission, though changes at lower energies (< 1 keV) were 
more important. While a higher density model would be more 
appropriate here, there is not yet one available for use. 

We use a model where the irradiating blackbody component 
is included in the reflection model, in order to calculate the 
reflection fraction, R, which is given in Table [8] with respect 
to that for a point source above a slab. The overall model 
is given by a combination of the reflected spectrum and this 
blackbody. In the model, the ionization parameter is defined 
as 



4tt£x 



(1) 



where, F x is the irradiating flux (erg cm s ) in the 0.001 
- 100 keV range and «h is the hydrogen number density. The 
ionization parameter, ^, therefore has units erg cm s . The 

2 The reflection models used here are to be made publically available on 
the XSPEC website 
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normalization of the reflection model is simply a scale factor 
that sets the overall strength of the combined blackbody and 
reflection spectrum, whereas the reflection fraction sets the 
relative strength of the two. The observed, unabsorbed flux 
(erg cm~ 2 s~ 1 ) of the blackbody plus reflection spectrum from 
the source in the 0.001 - 100 keV range is therefore given by 



^obs.unabs = Norm xF x (l+ R) 



(2) 



where, Norm is the normalization of the model. Rearranging 
gives 

4ft^obs,unabs 



Norm : 



ty H (l+R) 



(3) 



and the density used in the models, «h = 10 cm , is sub- 
stituted. 

To make these spectral fits more easily reproducible by 
other models we have calculated the equivalent blackbody 
normalization (using the definition of the XSPEC model 
bbody), from the best fitting value of the reflection model nor- 
malization and The strength of the reflection component 
relative to this is given by the reflection fraction. So that the 
results are reproducible using the same model as used here, 
we also quote the normalization value directly from the model 
(as defined above). 

The strength of the Fe K line is dependent on both the 
ionization parameter and the temperature of the irradiating 
blackbody. As the blackbody temperature increases, the num- 
ber of photons above the Fe K edge increases, leading to a 
greater EW. The EW peaks both at low ionization parameters, 
logc; = 1 .0, where Fe Koc is at 6.4 keV, and then at log^ ~ 2.7 
where the line is at 6.7 keV. At the peak for the higher ion- 
ization parameter, the EW is between ~ 100 - 300 eV for a 
blackbody temperature in the range 1 .5 keV to 3.0 keV. This is 
shown clearly in figure 3 of Ballantyne (2004). These equiva- 
lent widths are for a reflection fraction of 1 . Lower reflection 
fractions lead to smaller equivalent widths. 

We relativistically blur this reflection model by convolving 
with the rdblur model in XSPEC. rdblur is simply a con- 
volution with the diskline model, allowing the entire reflec- 
tion component to be blurred by the relativistic effects present 
at the inner accretion disk. This blurred reflection compo- 
nent is then included in the model instead of the diskline 
model used in the phenomenological fits. This modeling acts 
to self-consistently model the entire reflection spectrum (mul- 
tiple emission lines and continuum) rather than just treating 
the Fe K line alone. 

We are able to achieve good fits with this reflection mod- 
eling, generally finding similar inner disk radii as with phe- 
nomenological modeling (see spectral fit parameters in Tables 
[7] and [8). An example of our reflection fits is shown in Fig. [6] 
It is important to note that because the blackbody component 
drops off over 8-20 keV, so too does the reflection com- 
ponent. The typical strong Compton hump seen in the black 
hole sources (both in AGN and X-ray binaries) is therefore 
not apparent because of this. 

The two sources where we do not fit this black- 
body reflection model are the accreting millisecond pul- 
sars (SAX J1808.4-3658 and HETE J1900. 1-2455). Both 
sources show particularly hard spectra, where the high en- 
ergy continuum is dominated by a power-law. We therefore 
fit these spe ctra with the constant d ensity ionized disk reflec- 
tion model (Bal lantyne et alj |2001). where a power-law irra- 
diates a slab. As with the blackbody reflection model, the 
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FIG. 6. — Reflection model fit to the second Suzaku observation of 
GX 349+2 (left) and the second Suzaku observation of GX 17+2 (right). The 
individual model components are shown - red: disk blackbody, green: black- 
body, cyan: power-law, blue: reflection component. 

model normalization is dependent on and is given by the 
same definition as above. We therefore have calculated the 
equivalent power-law normalization at 1 keV so that these re- 
sults are reproducible using other reflection models. Again, 
we also report the best-fitting normalization directly from the 
model so that these results can be easily reproduced using 
the same model. This model fits the data well in both cases. 
HETE J1900. 1—2455 requires only the power-law reflection 
component, and no other continuum components are needed. 
SAX J1808.4-3658, on the other hand, requires both a disk 
blackbody and blackbody in addition to the power-law reflec- 
tion component. 

All spectral parameters from this reflection fitting are given 
in Tables |7]and[8] and source fluxes and luminosities are given 
in Table [6] We also show the range of measured inner disk 
radius from the reflection component in Fig. [4] We find a 
large number of source have their measured inner disk radii 
pegged at the model limit of 6 GMj c 2 , we discuss this in more 
detail later. 

5. SIMULTANEOUS TIMING OBSERVATIONS 

Fe K emission lines have been detected sim ultaneously 
with k Hz QPOs previously (e.g. GX 17+2, ICackett et at] 
2009b), however, there has not yet been the sensitivity to test 
whether the line and kHz QPO originate from the same region. 
While simultaneity is clearly important due to known changes 
in kHz QPO frequency with source state, comparisons of 
non-simultaneous measure ments of R m from F e K lines and 
kHz QPOs are promising (Cacket t et alJl2008h . Simultane- 
ous observations were performed by the RXTE Proportional 
Counter Array for both the Suzaku pointings of GX 17+2 and 
the second Cyg X-2 observations. Both sources are known 
kHz QPO sources. T hose QP Os were previous ly reported by 
IWiinands et all (fl99l andlHoman et all d2002b forGX 17+2, 
and by lWijnands et all dl998h for Cyg X-2. 

We have analyzed the PCA observations, considering seg- 
ments of continuous observation (ObsIDs): 2 and 5 ObsIDs 
worth of data were recorded for GX 17+2 and Cyg X-2 re- 
spectively. For each ObsID, we have computed an average 
Power Density Spectrum (PDS) with a 1 Hz resolution, using 
events recorded between 2 and 40 k eV. The PDS are normal- 
ized according to lLeahy et al.1 (119831) . so that in the absence of 
significant deadtime (as is the case here for both sources), the 
Poisson noise level is flat and constant with a value close to 2. 
The PDS so computed was then blindly searched for excess 
power between 300 Hz and 1400 Hz us ing a scanning tech- 
nique, as presented in Boi rin et all (2000). We failed to detect 
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TABLE 7 

Reflection model parameters: continuum 



Source 


Mission 


Obs. ID 


% 
10 22 cm" 2 


Serpens X- 1 


S k 


4U 1 U4ftU 1 u 


Ut/JT U.U_ 




AM!vl-:vt')\ lilt] 


inw inin in l 


U._1j it u.u 1 




XMM-New 'ton 


[ lUo-H nil 1 


U. jj it u.u 1 




XMM-Newton 


Any |mil(,n 1 
< II M- [\ l_U(H 1 J 


U._i / it U.U 1 


4U 1636—53 


XMM-New 'ton 


O.i 0_i J! 5 020 1 


0.37 ±0.01 




XMM-Newto n 


050035030 ] 


. 4 K ±0.01 




XMM-Ncwian 


0500350401 


0.47 ±0.01 


4U 1705-44 


Suzxiku 


401046010 


2.0 




Suzxiku 


401046020 


1.67 ±0.01 




Suzxiku 


401046030 


2.02 ±0.02 




XMM-Newton 


0402300201 


1.55 ±0.01 


4U 1820-30 


Suzxiku 


401047010 


0.21 ±0.01 


GX 17+2 


Stizxiku 


402050010 


2.10 ±0.06 




Suzxiku 


402050020 


2.15±0.05 


GX 349+2 


Suzxiku 


400003010 


0.77 ±0.01 




Suzxiku 


400003020 


0.80±0.01 




XMM-Newton 


0506110101 


0.67 ±0.01 


Cyg X-2 


Suzxiku 


403063010 


0.24 ±0.01 


SAX J1808.4-3658 


Suzxiku 


903003010 


0.30±0.01 




XMM-Newton 


0560180601 


0.27 ±0.02 


HETE J 1 900. 1 -2455 


Suzxiku 


402016010 


0.09±0.01 



Disk Blackbody 
Ar7- n (keV) 

1.22 ±0.01 
0.97 ±0.01 
1.01 ±0.01 
0.99±0.01 
0.21 ±0.01 
0.93 ±0.01 
1.00±0.01 
0.08 ±0.01 
1.04±0.01 
0.81 ±0.01 
1.42 ±0.01 
1.14±0.01 
1.69 ±0.02 
1.72 ±0.03 
1.55 ±0.01 
1.47 ±0.01 
1.24±0.01 
1.50±0.01 
0.41 ±0.01 
0.48 ±0.01 



Blackbody 



Norm. 


(keV) 


Norm (10 — 2 ) 




Norm. 


Line energy (keV) 


a (keV) 


Norm (10 — 2 ) 


110±5 


2.29 ±0.01 


5.3 ±0.5 


3.80±0.11 


1.19 ±0.01 








192 ±4 


1.79 ±0.01 


4.2±0.2 






1 .09 ±0.01 


0.13 ±0.01 


1.4 ±0.1 


166 ±6 


1.88 ±0.01 


2.95 ±0.07 






1 . 10 ±0.01 


0.17 ±0.01 


1.7 ±0. 1 


180± 1 


1.82 ±0.01 


3.53 ±0.06 






0.95 ±0.01 


0.29±0.01 


4.6±0.1 


1 0500 ± 2000 


1.92 ±0.02 


0.06 ±0.02 


1.94 ±0.06 


0.215 ±0.004 








78±6 


2.13±0.04 


0.9 ±0.3 


3.5±0.1 


0.51 ±0.01 








84±2 


2.10±0.01 


1.3 ±0.2 


3.5±0.1 


0.52±0.01 








..6+3-4) xlQ 8 


1.34 ±0.15 


00 , 6 +0.026 
UU -0.006 


1.68 ±0.01 


0.13±0.03 








107 ±4 


2.42 ±0.02 


2.3 ±0.4 












89±12 


2.12±0.01 


0.90 ±0.07 


3.02 ±0.04 


0.47 ±0.05 








2.63 ±0.1 1 


315-0.05 


0.22 ±0.06 












119 ±4 


2.50±0.04 


8.6 ±0.2 


2.38 ±0.05 


0.32 ±0.30 








99 ±5 


2.65 ±0.02 


6.7 ± 1.3 


1.8±0.5 


0017 +0.211 
1,01 -0.009 








95 ±4 


2.59±0.01 


9.4±1.3 


2.5±0.2 


27+ - 29 
U Z/ -0. 15 








109 ±4 


2.33±0.01 


11.5±2.1 












112±4 


2.35 ±0.01 


8.9±0.8 












221 ±6 


1.89 ±0.01 


15.1 ±0.1 






1.07±0.01 


0.09 ±0.01 


3.4±0.1 


163 ±1 


2.23 ±0.01 


3.5 ±0.2 


2.93 ±0.04 


0.87 ±0.09 


1.00 ±0.01 


0.11±0.01 


10 ± 1 


1370 ±120 


0.91 ±0.01 


0.05±0.01 


2. 14 ±0.02 


0.27±0.01 








232 ±21 


0.19 ±0.01 


0.6 ±0.1 


2.09 ±0.05 


0.31 ±0.01 














2.18±0.01 


0.07±0.01 









Note. — Continuum spectral filling parameters from reflection modeling. All uncertainties arc la. The x~ values for the fits are those in Tab.[s] Where no uncertainty is given, that parameter was fixed. Normalization of the disk blackbody is in 
fl 2 nner /(Z>10) 2 cosi where Rj n is the inner disk radius in km, /JlOis the distance to the source in units of 10 kpc, and i is the disk inclination. The normalization of the blackbody component is defined as L39/(D10) 2 , where L39 is the luminosity 

in units of 10 erg s — ^ . Normalization of the power-law component gives the photons keV - ' cm -2 s — ' at 1 keV. Finally, the normalization of the Gaussian is the total number of photons cm -2 s — ^ in the line. Note that the lower limit to the 
blackbody temperature in the reflection fits using the model ol Ballantyne 2004) (i.e. all except SAX J1808.4-3658 and HETE J1900.1— 2455) is 1.0 keV, and the upper limit is 3.15 keV. 



TABLE 8 
Reflection model parameters 



Source 




Mission 


Obs. ID 


Emissivily Index 


R iTl (GM/c 2 ) 


Blurred reflection parameleis 
Inclination (°) logE, 


Norm. (10 -26 


) Refl. fraction 


X 2 /dof 


Serpens X- 1 




Suzuka 


401048010 


2.2±0.1 


6+ 1 


16±1 


2.6 ±0.1 


13.6± 1.3 


0.18±0.01 


1526/1105 






XMM-Newton 


0084020401 


> 6 


107 ±11 


5±2 


2.8 ±0.1 


7.5 ±0.3 


0.10±0.01 


2213/1831 






XMM-Newton 


0084020501 


2.2 ±0.1 


15 ±2 


9±2 


2.6 ±0.1 


8.2±0.2 


0.16±0.01 


1982/1831 






XMM-Newton 


0084020601 


2.9 ±0.1 


44±2 


<3 


2.8±0.1 


6.4±0.1 


0.11 ±0.01 


2199/1831 


4U 1636-53 




XMM-Newton 


0303250201 


2.6 ±0.1 


8.5 ±2.0 


>68 


2.5±0.1 


0.20 ±0.05 


5-0-0.5 


1881/1832 






XMM-Newton 


0500350301 


3.2±0.4 


13±4 


63+j 6 


3.2 ±0.1 


0.54±0.20 


0.52 ±0.08 


2080/1832 






XMM-Newton 


0500350401 


2.9±0.1 


13±2 


>75 


3.0 ±0.1 


1.3±0.2 


0.42 ±0.02 


2339/1832 


4U 1705-44 




Suzxiku 


401046010 


3.0 


10.0 


24 


3.0±0.2 


00 , 7 +0.028 
UUi -0.007 


50-2.5 


1042/1127 






Sitzxikit 


401046020 


6.0 ±0.8 


7±1 


32 ±1 


3.0±0.1 


3.0±0.5 


0.43 ±0.02 


1387/1124 








401046030 


5.5 ±2.2 


11 ±6 


<4 


2.0±0.1 


8.7±0.7 


1.1 ±0.1 


1135/1122 






XMM-Newton 


0402300201 


>3.4 


34+f 2 


16±3 


3.0±0.1 


0.26 ±0.07 


0.21 ±0.02 


1766/1752 


4U 1820-30 




SltZuklt 


401047010 


>6 


11±1 


27±2 


2.6±0.1 


28.1 ±0.5 


0.11 ±0.02 


1336/1110 


GX 17+2 




Sitzxikit 


402050010 


2.2 ±0.1 


6+ 2 


< 11 


2.7±0.1 


15.5 ±2.9 


0.20 ±0.03 


647/606 






Siizuku 


402050020 


2.3 ±0.3 


6+ 2 


<6 


2.7±0.1 


20.0 ±2.7 


0.21 ±0.02 


684/606 


GX 349+2 




Sitztika 


400003010 


2.4±0.1 


10±2 


28±1 


2.9±0.1 


15.0±2.7 


0.24±0.01 


1341/1109 






SllZtlkl! 


400003020 


> 1.8 


386 ±97 


18±2 


3.1 ±0.1 


7.7±0.7 


0.22±0.01 


1276/1109 






XMM-Newton 


0506110101 


0.3 ±0.4 


6 +54 


48 ±2 


2.9±0.1 


22.3±0.1 


0. 14 ±0.01 


2500/1831 


Cyg X-2 




Suzuki! 


403063010 


2.0±0.1 


6+8 


<9 


2.6 ±0.1 


9.5 ±0.4 


0.31 ±0.03 


1379/1135 


SAX 11808.4 


-3658 


Sitztikit 


903003010 


>5.2 


14±2 


67_ 4 


2.7±0.1 


10.2±0.5 


0.86 ±0.08 


581/620 






XMM-Newton 


0560180601 


2.2±0.1 


6.0+ - 2 


67-1 


2.8±0.1 


8.6±0.2 


0.34 ±0.03 


2120/1952 


HETE J 1900. 


1-2455 


Snzakit 


402016010 


5.2±0.8 


14±1 


<4 


2.8±0.1 


2.1 ±0.3 


0.31 ±0.03 


1101/1227 


Note. — All uncertainties are lo. Where i 
Please see section[4*!2*ll'or further discussion. 


ro uncertainly i 


i given, lhat parameter ' 


>vas fixed. The n 


ormalizalion ol 


'the refleciion model (Norm.) scales 


ihe overall flux of thi 


> component. 



any kHz QPOs in any of the observations, for both sources. 
The analysis was also performed considering photons of en- 
ergy between 5 and 40 keV, but no kHz QPOs were detected. 
Finally, we have also produced an averaged PDS containing 
all observations. No QPOs were detected either. 

We have derived upper limits for QPOs of various widths, 
as listed in Table [9] for each pointed observation. Using data 
from combined observations, the upper limits are 15.5% and 
17.1% for GX 17+2, and 4.6% and 5.1% for Cyg X-2 for sig- 
nal of 100 and 150 Hz width respectively. These upper limits 
are not constraining for GX 1 7+2, whose lower and upper kHz 
QPOs have maximum RMS amplitudes of ~ 4% and ~ 6% 
respectively dWiinands et alj|1997t fHoman et alj|2002l) . This 
is explained by the combination of a very low source count 
rate and a reduced number of operating PCA units (2 or 3). 
For Cyg X-2, which displays larger source count rates, the 
RMS upper limits are significantly lower, but unfortunately 
still above th e maximum RMS amp litude of its kHz QPO re- 
ported so far (fWiinands et alJll998f) . 

Nevertheless for Cyg X-2, inspecting lower frequencies, 
features around 5 and 50 Hz are visible in the ObsID av- 



eraged PDS. These are likely normal branch and horizon- 
tal branch QPOs, previously known from the source. It has 
been shown in sources li ke Cyg X-2 (the p rototype being 
GX 17+2, Wij nands et alJfl997l Hoirian et al]p2002h that the 
upper kHz QPO frequency correlates with the HBO frequency 
up to a certain frequency where the upper QPO frequency 
keeps increasing while the HBO on e saturates (see for in- 
stance Figure 5 in lPsaltis etaT]|1999h . In Cyg X-2, the HBO 
frequency has been shown to increase from about 35 Hz up 
to a saturating frequency of 55 Hz. At the same time, the 
upper QPO frequency varied from ~ 730 Hz to ~ 1000 Hz, 
whi le the corresponding H BO frequency r anged from 35 to 55 
Hz dWiinands et al.| [l998). According to lPsaltis et aT] d 19991) 
(their fig. 8), before the saturation, a power-law relation- 
ship exists between the HBO and upper QPO frequencies, in 
the form of v H bo = 13.2ai(v upP er/1000) bl Hz, with a\ - 4.3 
and b\ - 1.4 in the case of Cyg X-2. The HBO QPO sig- 
nificance is ~ 5o significance (single trial). Fitting the 1 - 
2048 Hz PDS with two Lorentzians and a constant to account 
for the Poisson noise level, we have estimated the HBO fre- 
quency to be 49.5^i Hz, the QPO width as 7.0^ Hz for an 
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TABLE 9 

RXTE OBSERVATIONS OF GX 17+2, CYG X-2 AND HETE J 1 900. 1 -2455 



Name 




ObsID 


Date 


Time 


T bs 


Rate 


Bkg 


Npcu 


RMS 50 (%) 


RMS 100 (%) 


RMSiso (%) 


GX 17+2 




90022-07-06-01 


19/09/2007 


09:11 


2944.0 


79.5 


30.0 


2 


23.3 


28.7 


33.0 


GX 17+2 




90022-07-07-00 


27/09/2007 


15:15 


1392.0 


168.1 


53.7 


3 


17.5 


19.6 


21.3 


CygX-2 




93443-01-01-03 


01/07/2008 


01:34 


2016.0 


644.3 


40.8 


2 


6.8 


8.3 


9.4 


CygX-2 




93443-01-01-04 


01/07/2008 


03:12 


1776.0 


602.1 


37.0 


2 


5.8 


6.7 


7.6 


CygX-2 




93443-01-01-05 


01/07/2008 


04:50 


1536.0 


464.2 


39.3 


2 


8.1 


9.5 


10.5 


CygX-2 




93443-01-01-06 


01/07/2008 


06:28 


1408.0 


462.6 


37.4 


2 


6.6 


8.0 


9.0 


Cyg X-2 




93443-01-01-08 


01/07/2008 


11:09 


1344.0 


772.0 


16.5 


2 


6.0 


7.0 


7.9 


HETE J 1900.1 - 


-2455 


93030-01-16-00 


13/10/2007 


12:49 


3392.0 


88.3 


55.2 


3.0 


27.8 


34.9 


40.8 


HETEJ1900.1- 


-2455 


93030-01-17-00 


21/10/2007 


10:52 


3392.0 


120.3 


54.0 


3.0 


19.6 


24.1 


27.4 



NOTE. — T bs is the observation duration (s) and Rate and Bkg give the source and background count rates (c/s) respectively. The number of active 
PCA units (N PCU ) are listed for each ObsID. The upper limits on the QPO RMS amplitude (%) have been computed for QPO width of 50, 100 and 150 Hz 
(RMS50, RMS100, RMS150 respectively). They have been computed at the 3a confidence level, meaning that the RMS listed would enabl e the detection 
of any QPOs with a 3a accuracy on its amplitude (this corresponds to a ~ 6a excess (single trial) in the power density spectrum, see Boutelier et al. 
1 2009) for details). 



RMS amplitude of about 3.5^ 5%- The latter two parameters 
are fu lly consistent with previous detections dWiinands et alj 
1998). Despite its low significance, it is therefore tempting 
to use the ~ 50 Hz HBO reported here (which is below the 
saturating frequency) to estimate with the above formula, the 
frequency of the upper QPO, that we failed to detect. This 
yields a frequency of about 910 Hz. 

As we noted earlier, if both the upper kHz QPO and Fe K 
emission line arise from the same region in the disk, and the 
kHz QPO is associated with an orbital freque ncy, they can 
be combined to estimate t he neutron star mass dPiraino et alj 
2000; ICackett et alJl2008h . If we tentatively use the inferred 
910 Hz from scaling of the HBO frequency, along with the 
inner disk radius from the iron line we get a mass estimate 
in the range 1.5-2.4 M (depending on whether we use the 
phenomenological or reflection mo deling results), consistent 
with the optical ma ss measurement (Orosz & Kuulkersl l 19991 : 
Elebert et al. 2009). This implies a consistent radius between 
the Fe K emission line and the kHz QPO. But, again, we are 
extremely cautious about this given the large uncertainty in 
scaling the QPO frequency, and the possible range of inner 
disk radius. 

Looking at the RXTE archive, two observations of 
HETE J 1900. 1—2455 were performed contemporaneously 
with the Suzaku pointings of this source. Unfortunately, no 
kHz QPOs were detected. Upper limits are also listed in Ta- 
ble [9] but they are not very constraining on the presence of 
kHz QPOs due to the low source count rate and a compara- 
tively large background. 

Archival RXTE observations are also available for the three 
XMM-Newton p ointings of4U 1636 —53: the first observation 
was analyzed in Bar ret et all (12 005) and the latter two are also 
considered in lBoutelier et alj (2010). We analyze the data us- 
ing the same method as described above. As strong kHz QPOs 
were detected in many pointed observations, the QPO param- 
eters (mainly the quality factor) we re corrected for freque ncy 
drift using the method described in lBoutelier et alJ (120091) . A 
total of 10 observations were considered. In 4 of those, we 
detect twin kHz QPOs, in one we detect a single upper kHz 
QPO at — 500 Hz (i.e. the lower kHz QPO which is expected 
at much lower frequencies has merged with the noise and is 
no longer a QPO), and in two, a single lower kHz QPO (for 
an identification of single kHz QPO, b ased on th eir quality 
factor and RMS amplitude see, iBarret et all 120061) . The non 
detection of the upper kHz QPO in the latter t wo pointed ob- 
servati ons is due to a lack of sensitivity, see iBoutelier et alj 
(2010), but the frequency of the upper kHz QPO can be ap- 



proximated by t he nearly linear fu nction that links the two 
frequencies (e.g. B elloni et al.ll2005l) . The results of the fits 
are listed in Table [TOl 

These simultaneous XM M-Newton and RX TE observations 
are discussed extensively in Altam irano et alJ (120101) : here we 
briefly discuss the consis tency of our spectral fitti ng with the 
kHz QPO properties. As Altamira no et alj ((2010) find, when 
assuming a 1 .4 M neutron star, the inner disk radius (as mea- 
sured by the Fe K line from phenomenological fits) disagrees 
with the implied inner disk radius assuming that the upper 
kHz QPO is a Keplerian frequency in the disk. If instead one 
uses the phenomenological Rj n measurements combined with 
the upper kHz QPO frequency to estimate a mass, the implied 
mass would be 4.2 ±0.5 M , 1.7 ±0.2 M and 2.2 ±0.1 
M from the three observations respectively (we use the kHz 
QPO frequencies from the RXTE observations that overlap 
with the XMM-Newton observations for the longest period). 
Clearly, the first observation estimates far too high a neutron 
star mass, while the other two observations give a more rea- 
sonable (though quite high) estimate. Note a high mass is 
predic ted for 4U 1636— 53 from the drop in kHz QPO coher- 
ence (IBarret et al.ll2.006T> . From reflection modeling, we find 
slightly different R m measurements than for the phenomeno- 
logical fits. The corresponding mass estimates from reflection 
fits are: 2.6 ± 1.0 M , O.8±O.4M andO.7±O.2M . Again, 
the separate measurements are not all consistent with each 
other, and this time two of the masses are clearly on the low 
side. 

Looking at the trends, if the two phenomena originate in 
the same part of the disk, a lower kHz QPO frequency should 
correspond to a larger inner disk radius (as measured by the 
Fe K line). However, in 4U 1636—53 the opposite behavior 
is seen - as the kHz QPO frequency increases there is not 
a corresponding decrease in R m , in fact, when using the R m 
measurements from reflection modeling R m is actually seen to 
increase. This, however, is only based on three observations, 
and further, higher quality, observations are needed to make 
any strong conclusions. 

6. DISCUSSION 

We have presented a comprehensive study of Fe K emission 
lines in 10 neutron star LMXBs. In all cases, the emission line 
profiles are well modeled by a relativistic disk line assuming 
a Schwarzschild metric. However, in one source, GX 340+0, 
the exact line profile is uncertain: a relativistic emission line 
fits the data well but only if either the abundance of Fe is non- 
standard, or alternatively there is an Fe XXVI absorption line 
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TABLE 10 

RXTE OBSERVATIONS OF 4U 1636-536 SIMULTANEOUS WITH THE THREE XMM-Newton OBSERVATIONS 



Name 



ObsID 



Date 



Time 



1 obs 



Rate Bkg N P , 



RMS (%) 



4U 1636-53 
4U 1636-53 
4U 1636-53 
4U 1636-53 



70034-01-01-01 
91027-01-01-000 
91027-01-01-000 
93091-01-01-000 



29/08/2005 03:24 

29/08/2005 16:35 

29/08/2005 18:06 

28/09/2007 14:47 



4338.0 
4434.0 
23362.0 
4464.0 



406.9 53.3 

420.6 78.8 
347.4 78.9 

482.7 40.6 



4U 1636-53 93091-01-01-000 28/09/2007 16:17 23392.0 415.7 36.8 



4U 1636-53 
4U 1636-53 
4U 1636-53 

4U 1636-53 

4U 1636-53 



93091-01-01-00 
93091-01-02-000 
93091-01-02-000 

93091-01-02-000 

93091-01-02-00 



28/09/2007 22:47 
27/02/2008 03:46 
27/02/2008 04:35 



7104.0 
1280.0 
4352.0 



179.5 17.1 
690.2 58.3 

688.6 34.9 



27/02/2008 06:07 20352.0 546.4 31.0 
27/02/2008 11:46 10096.0 252.5 17.0 



509.6 ±35.6 



643.5 ±2.5 
984.1 ±19.2 
589.7 ±10.5 
906.3 ±5.9 

645.9 ±5.7 
646.0 ±2.6 
979.5 ±27.1 
701.6±3.7 
1018.8±32.8 
882.0 ±2.6 



191.4±85.8 12.3±2.0 3.0 



17.8±8.3 
165.3±67.8 
94.6±31.8 
88. 9 ±15.8 

36.5 ±15.8 
19.3±7.4 
184.9±71.8 
86.2±13.5 
221.3± 111.4 
15.8±5.8 



6.5 ±1.0 3.4 

U.9±1.6 3.8 

7.8±0.9 4.2 

10.0±0.6 7.7 



8.4±1.2 
5.5 ±0.7 
9.3 ± 1.3 
10.5 ±0.6 
7.8± 1.4 
6.7±0.9 



3.4 
3.8 
3.6 
8.5 
2.8 
3.7 



NOTE. — T b* is the observation duration (s) and Rate and Bkg give the source and background count rates respectively. The number of active PCA 
units (Npcu) are listed for each ObsID. The frequency, Full Width Half Maximum, and RMS amplitude (%) of each kHz QPO are reported with their lo 
error. R is the ratio of the Lorentzian amplitude to its lo" error. An R of 3 corresponds roughly to 6o" excess power in the Fourier PDS, see Boutelier et al. 
(2009) for details. When two QPOs are detected in one ObsID, the first one refers to the lower kHz QPO, whereas the second one refers to the upper kHz 
QPO. 



superposed on the emission line. Fitting of broadband spec- 
tra, makes it clear that the hot blackbody component domi- 
nates the spectrum from around 7-20 keV (in the Z and atoll 
sources), and thus provides the vast majority of the flux irra- 
diating the accretion disk and leading to the Fe K emission 
line. We fit reflection models to the spectra, finding that in the 
Z and atoll sources a reflection model where the irradiating 
source is a blackbody does fit the spectra well. In the ac- 
creting millisecond pulsars, however, a power-law dominates 
above 10 keV, and standard reflection models used to fit black 
hole spectra work well there. 

6.1. Inner disk radii 

The inner disk radii that we measure mostly fall into a small 
range: 6-15 GM/c 2 . Also consistent with this range is 
a recent independent measure ment of the inner di sk radius 
in 4U 1705-44 (R in = 14 ± 2: IPiSaivo et alJl2009h . Where 
there are multiple observations of a source, in most cases there 
is not a large change in inner disk radius, which may suggest 
the inner disk radius does not vary with state. While we note 
this from observations at different epochs, similar results have 
been found previously when looking at how the line profile 
changes with source state during one observation (iD'Ai et alJ 
2009: IIariaet aT]|2009h . 

The only source where we see variations in inner disk radii 
is Ser X-l, and this is easily apparent when one compares 
the line profile observed by Suzaku with one of the XMM- 
Newton observations (see Fig. |3j. Comparing these observa- 
tions, the spectrum from the Suzaku observation of Ser X- 
1 is harder than the XMM-Newton observations - although 
both the blackbody and disk blackbody components are hot- 
ter, the relative normalization of the blackbody component 
relative to the disk blackbody component is higher. The ion- 
izing flux during the Suzaku observation is therefore stronger, 
which would explain the higher line flux. In all cases, though, 
the emitting radius inferred from the blackbody component 
is very similar. The 0.5 - 10 keV source flux is also high 
est during the Suzaku observation (10.3 x 10 
compared to 6.3 x 10~ 9 ,4.4 x 10 

s _1 for the XMM-Newton observations). This may tentatively 
suggest that the disk extends closer to the stellar surface at 
the highest luminosities, though this is not confirmed when 
looking at all sources (Fig. [7). 



erg cm 2 s 1 
9 and 4.9 x 10~ 9 erg cm~ 2 



™o30 



~20 



uio - 




Luminosity (erg/s) 

FIG. 7. — Measured inner disk radii from phenomenological fits as a func- 
tion of 0.5-25 keV source luminosity. The dotted line indicates 6 GM/c 2 , the 
smallest allowed radius in the model. 



We have also looked to see if there is any significant change 
in inner disk radius with 0.5 - 25 keV luminosity for all 
sources (see Fig. [7). There is no apparent systematic trend 
with luminosity, and the inner disk radius is found to be close 
to 6 GM/c 2 for the majority of sources/observations over al- 
most 3 orders of magnitude in luminosity. Similarly, we see 
no clear difference between atolls, Z sources and accreting 
millisecond pulsars, as can be seen in Fig. |4] 

There is typically consistency between the inner disk radii 
measured from phenomenological and reflection spectral fits. 
Nevertheless, in both cases many of the inner disk radii 
measurements are pegged at the lower limit of the model 
(the innermost s table circular orbit, ISCO). In 4U 1705-44, 
iReis et aD (l2009bl) found that fitting with a model for a Kerr 
metri c (relevant for maximally spinning black holes lLaorl 
[1991 gave inner disk radii much smaller than expected neu- 
tron star radii. They suggested that as the disks in these sys- 
tems are ionized one must also take broadening due to Comp- 
tonization into account by using relativistic ionized reflec- 
tion models. When they did this, the inner radius measured 
was significantly higher, and no longer smaller than expected 
stellar radii. The reflection models employed here also take 
Compton broadening into account, yet, we still find that the 
inner disk radius pegs at the ISCO in a number of sources. 

This may not necessarily be a problem. As neutron stars 
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are spinning, their angular momentum will change the radius 
of the ISCO m aking it slight ly smaller for a corotating disk 
(see fig 1(a) of Miller 2007, to see how the ISCO changes 
with angular momentum). Potentially, one could then mea- 
sure the dimensionless angular momentum parameter, under 
the assumption that the inner disk is truncated at the ISCO 
and not at the boundary layer or stellar surface. To test this, 
we fit a line profile with a variable an gular momentum pa- 
rameter dBrenneman & Re ynolds 2006) to several of the best 
datasets in our sample. However, in each case we found that 
this parameter was not well constrained, being consistent with 
zero. Better, more sensitive data may be able to constrain this 
further. 

We can also compare the inner disk radii from the Fe line 
fitting with that implied from the disk blackbody fits. The nor- 
malization for the disk blackbody used here (diskbb) is de- 
fined asN = /? 2 nner cos //(D10) 2 , where ^i nner is the inner disk 
radius in km, D10 is the source distance in units of 10 kpc and 
i the disk inclination. Therefore, from the fitted normalization 
and the inclination determined from the Fe line fits, we can 
calculate the inner disk radius in km. Assuming a reasonable 
range for the neutron star mass (here we use 1.4-2.0 M Q ) we 
can convert to units of GM/c 2 for direct comparison with the 
Fe line fits. The disk radii determined this way are given in 
Table fTTI Where only an lower/upper limit for the inclination 
is found from the Fe line fits, we just use this limit for the 
inclination in the calculation. 

Comparing the inner disk radii from the continuum 
(diskbb) and Fe line fitting, we find that the the continuum 
inner disk radius is generally smaller than the Fe line value, 
though there are a small number of cases where the values are 
consistent (when assuming a mass of 1.4 M Q ). It is important 
to note that we have not used any spectral hardness factor or a 
correction for the inner boundary condition for the continuum 
radii. These two factors are expected to lead to a factor of 
1.18 - 1.64 increase between the ra dius measured by d iskbb 
and the real inner disk radius (e.g. Kub ota et al.l l2001). This 
would bring more of the values to be consistent, though a large 
number would still disagree. Where there are large disagree- 
ments, it could be possible that the continuum model chosen 
is not physical. For disagreements where the continuum ra- 
dius is smaller than the Fe line radius, one can imagine that 
disk photons may be scattered out of the line of sight by a 
corona, diminishing the strength of the corona, and leading to 
a smaller measured continuum radius. Additionally, the con- 
tinuum method relies on knowing the absolute flux of the disk 
component. Therefore, if the absolute flux calibration of the 
telescope, or if the continuum model (column density or other 
parameters) are incorrect, then this can lead to an incorrect 
inner disk radius. 

6.2. Source inclination and disk emissivity 

We found in many cases that there are several local % 2 min- 
ima within the diskline parameter space. There is usually 
one solution that is clearly the global minimum, with a sig- 
nificantly better % 2 value. However, we found that in one in- 
stance (GX 349+2, second Suzaku observation) the difference 
was as small as A% 2 ~ 2. This issue is due to a degeneracy 
where a model with a combination of low inclination, high 
emissivity index and high line energy can look similar to a 
higher inclination, lower emissivity index, and lower line en- 
ergy model. Usually the inner disk radius remains consistent 
between the solutions. This is apparent in Tab. [5] see partic- 
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FIG. 8. — Dependence of the x 2 confidence contours on inclination and 
emissivity index for GX 349+2. The cross marks the global x 2 minimum 
when the second Suzaku observation is fit on its own, and the star marks the 
minimum when all three observations are GX 349+2 are fit jointly, with the 
inclination tied between them. The contours are the lo~ (black), 90% (red), 
and 99% (green) confidence levels, with the solid lines from fitting the second 
observation on its own, and the dotted lines from when all three observations 
are fit jointly. 

ularly the GX 349+2 and Ser X-l observations. We demon- 
strate this issue by showing the dependence of the % 2 con- 
fidence contours on inclination and emissivity index for the 
second GX 349+2 Suzaku observation in Fig. [S] Here, one 
clearly sees two minima with only a small difference in % 2 . 
The % 2 minimum is 1262.6 for the low emissivity index, high 
inclination minimum compared to 1265.0 for the high emis- 
sivity index, low inclination local minimum (both have 1106 
degrees of freedom). For the first Suzaku GX 349+2 we also 
find two minima. In this case, however, the global minimum 
has high emissivity index, low inclination opposite to the sec- 
ond observation (% 2 = 1268.1 for low inclination compared to 
1273.6 for higher inclination). 

As discussed, comparing all observations of GX 349+2, the 
inclination and emissivity index from spectral fitting varies 
significantly (see Tab. [5]). This was first pointed out by 
llaria et aITj 2009) when compa ring their XMM-Newton results 
with the ICackett et"ai1 (120081) results. As they also noted, 
while it is possible that emissivity index and inner disk radius 
and the line energy may vary with source state, the source in- 
clination obviously does not. As we have demonstrated, this 
difference is not because the lines look very different. In fact, 
if all three observations of GX 349+2 are fit simultaneously 
with the inclination tied between observations (best fit value 
i = 44^4 degrees), then a good fit is achieved, and the lines 
look almost identical (see Fig. [9]). Plainly, a priori knowledge 
of source inclination is of significant benefit. However, joint 
fitting of multiple observations of the same source potentially 
resolves this issue. As can been seen in Fig. [8] when the three 
observations of GX 349+2 are fit jointly, there is no longer a 
local % 2 minimum with a high emissivity index and low incli- 
nation. 

Such a degeneracy between the line parameters is self- 
evident when studying how each parame ter affects the line 
profile (see fig. 1 of Fa bian et al.l I1989I) . The inclination 
shifts the peak of the blue wing of the line profile. There- 
fore, a higher line energy can be offset by a lower inclina- 
tion. Moreover, a higher inclination makes the line broader 
and less 'peaked', and the emissivity index also has a similar 
effect. Thus, a higher inclination, and lower emissivity index 
can look similar to a lower inclination and higher emissivity 
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TABLE 1 1 
Disk blackbody inner radius 



Source 


Mission 


Obs. ID 




Phenomenological 




Reflection 










Staler (km) 


R in (GM/c 2 ) 
for 1.4 M Q 


R in (GM/c 2 ) 
for 2 M e 


Rinncr (km) 


R in (GM/c 2 ) 
for 1.4 M 


R m (GM/c 2 ) 
for 2 M 


Seipens X- 1 


Suzaku 


401048010 


8.7 ±2.2 


4 2 zt 1 


2 9 ± 7 


9.0±2.3 


4.4± 1.1 


3.0±0.8 




XMM-Newton 


0084020401 


13.1 ±3.2 


6 3 zb 1 6 


4 4 ± 1 1 


1 1.7 ± 2.9 


5.6±1.4 


4.0 ±1.0 




XMM-Newton 


0084020501 


U.3±2.8 


5 5 zb 1 4 


3 8 ± 1 


10.9 ±2.7 


5.3± 1.3 


3.7 ±0.9 




XMM-Newton 


0084020601 


12.7 ±3 J 


6.1 i 1.5 


A 1 -I- 1 1 
4. 5 ±1.1 


11.3 ±2.8 


5.5± 1.4 


3.8 ± 1.0 


4U 1636-53 


XMM-Newton 


0303250201 


143 ±72 


70 zb 35 


49 ± 25 


100 ±49 


49 ±24 


34 ± 17 




XMM-Newton 


0500350301 


12.7 ±6.3 


6 2 zb 3 


4 3 ± 2 1 


7.9±2.2 


3.8 ±1.1 


2.7 ±0.7 




XMM-Newton 


0500350401 


18.8±9.4 


9.1±4.6 


6.4±3.2 


10.8 ±5.3 


5.2±2.6 


3.7 ± 1.8 


4U 1705-44 


Suzaku 


401046010 


291(1+4560 


1400+ 2200 
i+uu_ 490 


gon+1540 


9780+ 6410 


4740+ 3100 
+ '+u_ 174 o 


3320+ 21 ™ 




Suzaku 


4U 1 U40UZU 


+ ..> ± U.Z 


2.2±0.1 


1.5±0.1 


O.J ± U. J 


J.Z ± U. I 


Z.Z ± U. 1 




Suzaku 


401046030 


5.3 ±0.2 


2.6 ±0.1 


1.8 ±0. 1 


5.5 ±0.4 


2.7 ±0.2 


1.9 ±0. 1 




XMM-Newton 


0402300201 


1.7 ±0.9 


0.8 ±0.4 


0.6 ±0.3 


0.96 ±0.04 


0.47 ±0.02 


0.33 ±0.01 


4U 1820-30 


Suzaku 


401047010 


7.8±0.5 


3.8±0.2 


2.6 ±0.2 


8.8±0.5 


4.3±0.2 


3.0±0.2 


GX 17+2 


Suzaku 


402050010 


9.2 ±0.6 


4.5 ±0.3 


3.1 ±0.2 


9.8 ±0.5 


4.8 ±0.2 


3.3±0.2 




Suzaku 


402050020 


10.0±0.5 


4.8 ±0.2 


3.4±0.2 


9.6 ±0.4 


4.6 ±0.2 


3.2±0.2 


GX 349+2 


Suzaku 


400003010 


5.6±1.4 


2.7 ±0.7 


1.9 ±0.5 


5.6±1.4 


2.7 ±0.7 


1.9 ±0.5 




Suzaku 


400003020 


5.9± 1.5 


2.8 ±0.7 


2.0 ±0.5 


5.4±1.4 


2.6 ±0.7 


1.8 ±0.5 




XMM-Newton 


0506110101 


10.6 ±2.7 


5.1 ±1.3 


3.6 ±0.9 


9.1 ±2.3 


4.4± l.l 


3.1 ±0.8 


Cyg X-2 


Suzaku 


403063010 


14.8±2.7 


7.2± 1.3 


5.0±0.9 


14.1 ±2.6 


6.8±l.2 


4.8 ±0.9 


SAX J 1808 


Suzaku 


903003010 


1.1 ±0.1 


0.51 ±0.04 


0.36 ±0.03 


20.7 ±2.0 


10.0± l.O 


7.0±0.7 




XMM-Newton 


0560180601 


ntf 


34^' 


24+» 


8.5 ±0.5 


4.l±0.2 


2.9 ±0.2 



NOTE. — Assumed distances are given in Table [Tl 



high emissivity index, causing the irradiation of the disk to be 
more centrally concentrated. In the most extre me case of spin- 
ning black holes an index of ~ 5 is expected (Miniutti et al. 
2003). These effects should not be as strong in neutron stars, 
and thus a lower index would be expected. 

It is also important to note that in the reflection modeling 
the inclination is only introduced by the relativistic blurring, 
and the reflection spectrum itself is angle-averaged. However, 
reflection is sensitive to angle, with the reflection strength 
weakening with higher inclination. The use of an angle- 
averaged model may therefore bias the results obtained. 

In black hole systems, the ionizing flux is dominated by a 
power-law. The Compton back-scattering part of the reflec- 
tion spectrum is easily detected against a power-law contin- 
uum. Disk reflection is more difficult to detect against a ther- 
mal spectrum since the continuum flux drops precipitously 
with higher energy. In the case of neutron stars, then, disk 
reflection modeling may only yield improved physical con- 
straints in the limit of extraordinary high energy sensitivity. 
The fact that many of our reflection fits imply small radii 
and low inclinations may be partially due to a combination 
of modest high energy sensitivity, coupled with the limita- 
tions of angle-averaged disk reflection models. Deeper obser- 
vations with Suzaku aimed at clearly detecting the Compton 
back-scattering hump will help to obtain better physical con- 
straints. 

6.3. Ionization parameter and reflection fraction 

Reflection modeling allows us to infer the ionization state 
of the inner accretion disk. The reflection spectrum is highly 
dependent on th e ionization parameter, Z, (see fig. 1 of 
Ballantyne 2004). For low ionization states many emission 
lines are present which disappear as the ions become more 
highly ionized. The Fe Koc line in particular, evolves signif- 
icantly with ionization parameter. For neutral reflection, the 
6.4 keV line is strong and as the ionization increases this line 
becomes weaker. The 6.7 keV line then becomes prominent 
until it disappears as iron becomes fully ionized. This evo - 
lution is displayed clearly in figure 3 of lBallantvnd ([2004), 
where the dependence of the equivalent width (EW) of the 



T 



1.1 - 




Energy (keV) 

FIG. 9. — Fe K line profile (shown as data/continuum model) for the three 
observations of GX 349+2. Here the inclination was tied between the obser- 
vations during the spectral fitting. Black shows the first Suzaku observation, 
green the second, and red the only XMM-Newton observation. 

index. 

We find a similar problem with the inclination determined 
from fitting the Cyg X-2 spectra, where we measure an in- 
clination of 25°. This inclination is not consistent with the 
value measured from optical observations, where lower and 
upper inclinatio n limits are 49° - 73 ° , with a best-fittin g 
value of 63° dOrosz & Kuulkerslll999t lElebert et al.ll2009l) . 
If we constrain the inclination to be within the optical val- 
ues we do not achieve as good a fit (% 2 — 1337 compared to 
% 2 = 1305), and find the inclination pegs at the lower bound 
of 49°. With the inclination constrained, the line parameters 
we find are E = 6.40+ 01 keV, q = 2.9 ± 0.2, R in = 25 ± 4 
GM/c 2 , i — 49 +1 °. Again, there is a pattern where a higher 
inclination leads to a lower line energy and emissivity index. 

Note that in several cases we find high emissivity indices 
and/or very low inclinations when fitting the reflection mod- 
els. Both a high emissivity index and low inclination try to 
make to line profile more peaked. We tested whether this can 
be solved with a higher Fe abundance by using a reflection 
model with two times solar abundance, but retrieved similar 
results, and worse fits. Light bending can potentially lead to a 
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Fe Ka on ionization parameter is shown. The EW peaks for 
either neutral reflection when the 6.4 keV line strong, or for 
log^ = 2.6 - 2.8 when the 6.7 keV line is dominant. From 
our reflection fitting (Tables [7] and [8) we find quite a narrow 
range of ionization parameter. Most spectra are best fit with 
an ionization parameter in the range logt, = 2.6 - 2.8, where 
the Fe Ka EW peaks. 

We also consider the strength of the reflection component, 
which directly relates to the solid angle subtended by the re- 
flector, CI. The reflection fraction, R — H/2%, is given in Table 
[8] A reflection fraction of 1 corresponds an isotropic source 
above the slab. Most of the observations here have a reflection 
fraction in the range R = 0.1 - 0.3, suggesting that a reason- 
able fraction of the ionizing flux is intercepted by the accre- 
tion disk and re-emitted. If we interpret the blackbody com- 
ponent as the boundary layer, this may also suggest that the 
boundary layer has quite a similar extent in th ese sources. As 
discussed above, figure 3 of Ballantyne (2004) shows a strong 
dependence of the Fe Ka line EW and the blackbody tempera- 
ture. Most of the blackbody temperatures we find have kT > 2 
keV, which equates to an EW of > 200 e V for R = 1 . From the 
phenomenological fits (see Table|5j, we get equivalent widths 
that are mostly less than 100 eV, thus, the reflection fractions 
that we find generally seem reasonable. 

There are two cases (one spectra of 4U 1705—44 and one 
spectra of 4U 1636—53) where the best-fitting reflection frac- 
tions are at the upper limit of the model (R = 5). In these 
spectra, either the blackbody temperature or the blackbody 
normalization is very low. Both these lead to a low equiva- 
lent width, and therefore for the model to achieve the required 
line strength a large reflection fraction is required. Both these 
sources a re atolls which do trans it between hard and soft states 
(see, e.g-. lGladstone et aT1l2007l) . and these two particular ob- 
servations are in the hard state. As discussed earlier, spec- 
tral deconvolut ion in atoll and Z sources is not unique (e.g. 
iLin et al.l 120071) . and the high values of the reflection frac- 
tion suggest that our particular choice of model here is un- 
likely the best physical model. We have also tried replacing 
the blackbody reflection model with a power-law reflection 
model. While this achieved more reasonable refl ection frac- 
tions i t did not provide as good a fit to the data. ID' Ai et aD 
(2010) have recently successfully modelled the hard state in 
4U 1705-44 using an alternative model involving a Comp- 
tonization model to fit the hard component, and such a model 
may be more appropriate for these two particular observa- 
tions, though we do not investigate this here. 

6.4. The boundary layer illuminating the disk 

It is well known that a hot blackbody (or alterna- 
tively a Comptonize d componen t ) dom inates from around 
7-20 keV (e.g. iBarret et al.l 120001) . In all sources 
but the accreting millisecond pulsars, our spectral fitting 
agrees with this. This blackbody component may po- 
tentially be due to emission from the boundary layer be- 
tween the inner accretion disk and neutron star surface 
(see, e.g.. ISunvaev & Shakuralll986l; llnogamov & Sunvaevl 
1999HPopham & Sunvaevll200UlGrebenev & Sunvaevll2002t 
Gilfanov et al.l2003tlRevnivtsev & Gilfanovl2006l fordiscus- 
sions about the boundary layer). The boundary layer is the 
transition region between the rapidly spinning accretion disk 
and the more slowly spinning neutron star, and should account 
for a significant fraction of the accretion luminosity. 

As neutron star spectra show that the blackbody compo- 
nent dominates the flux able to ionize iron, it suggests that the 



boundary layer may be the dominant source of ionizing flux 
illuminating the inner accretion disk, leading to Fe Ka emis- 
sion. Our spectral analysis with reflection models shows that 
this picture is consistent with the data. 

It is not a new idea that the boundary layer or neutron star 
surface c an illuminate the dis k leading to Fe K emission. For 
instance, lDav&Donel(fT99Tl) were the first to note that X-ray 
bursts on the surface of the neutron s tar can lead to a disk- 
reflection component. Furthermore, Bran dt~& Matt! (119941) 
use a geometry where an inner spherical cloud of Comp- 
tonized gas acts to illuminate the inner accretion disk as a 
basis for calculating Fe K profiles expected in Z sources. 
Moreover, d etailed calculations of the in ner disk and bound- 
ary layer led Popha m~& Sunvaevl (120011) to conclude that the 
X-ray flux from the hot boundary layer incident on the accre- 
tion disk would be large enough to lead to line emission and 
a reflection spectrum. 

Observationally, there is not a large body of previous work 
fitting broadband neutron star spectra with reflection models 
wher e the boundary lay er irradiates the accretion disk. How- 
ever, |DoneeraL| (|2002|) successfully fit Ginga data of Cyg X- 
2 with a reflection model which illumi nates the disk with 
a Comptonized continuum, and similarly Gierlihski & Pond 
( 120021) achieved good fits to the RXTE spectra of atoll 
4U 1608—52 with the same model. However, the data in both 
cases had sig nificantly lowe r spectra l resolution to those pre- 
sented here. iDi Salvo et al.l d2000al) also successfully fit the 
broadband BeppoSAX spectrum of MXB 1728-34 with a re- 
flection model which illu minates the d i sk wit h a Comptonized 
continuum. Furthermore, Barr et et al.l yOOO) suggest that the 
boundary layer provides the X-ray flux irradiating the disk 
from s pectral fitting of RXTE data. 

The Popham & S unvaevl (120011) model for the boundary 
layer is one where t he hot boundary laye r gas is radially and 
vertically extended. Popham & SunvaevJ (120011) showed that 
the hot boundary layer is much thicker than the inner disk 
(see their fig. 1). As mass accretion rate changes they find 
that both the radial and vertical extent of the boundary layer 
will change, and thus the fraction of X-ray emission from the 
boundary layer that is intercepted by the disk will vary. An 
alternative model for the boundary layer is one where mat- 
ter spreads over the neutron star surface in a latitudinal belt 
(Inogamov & Sunvaevl ll999l) . In this spreading layer the belt 
width is dependent on the mass accretion rate. At low mass 
accretion rates (3 x 10~ 3 ) the width is of the order of the thin- 
disk thickness, however, the belt disappears at high Eddington 
fractions (~0.9) as the entire surface is radiating. It therefore 
seems reasonable that even the spreading layer may plausibly 
illuminate the inner accretion disk. 

We now further investigate the notion that the boundary 
layer is the source of ionizing flux. Using ionization parame- 
ters typical to what we find from spectral fitting, we determine 
the maximum height of the boundary layer, for a disk extend- 
ing close to the neutron star surface. The ionization parameter 
at the inner accretion disk, is given by 



where Lbl is the boundary layer luminosity, n is the number 
density in the disk, and R is the distance from the ionizing 
source to the disk. For the purposes of this estimation we use 
the distance between the top of the boundary layer and the 
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inner accretion disk. That distance is given by 

R = V^in+Z 2 ( 5 ) 

where R m is the inner disk radius, and Z is the height of the 
ionizing source above the disk. The distance will, of course, 
be higher for other regions of the accretion disk, and for 
boundary layer emission from closer to the equator. However, 
combining these two equations and re-arranging leads to the 
height of the ionizing source above the disk being defined as 



From spectral fitting and reflection modeling we determine 
Lbl, Rm an d an d from standar d disk theory we can put rea - 
sonable constraints on n. From Shakura & Sun vaevl d 19731) . 
we find n > 10 21 cm~ 3 . Forlog^ = 3, Lbl = 10 37 erg s _1 and 
Rya = 20 km, we find Z < 24 km. It is therefore completely 
reasonable that the boundary layer provides the ionizing flux 
to produce the Fe K emission lines observed. A subsequent 
detailed analysis of multiple spectra from 4U 1705—44 has 
also concluded that the bounda ry layer likely illu minates the 
disk, at least in high flux states dD'Ai et al. 2010). 

A further check can be performed by looking at the inferred 
blackbody emitting radius from the spectral fits. In Table [T2l 
we show the calculated emitting radius and assumed distances 
from the blackbody normalizations. We use the same dis- 
tances as given in Table [TJ As the boundary layer emission is 
more correctly described by a Comptonized component as op- 
posed to a simple blackbody we can maybe not interpret these 
numbers directly. However, they do provide a good consis- 
tency check, with radii typically a signi ficant fraction of the 
stella r radius, as expected theoretically ( Popham & Sunvaevl 
lIOOll) . In a couple of cases the radii are extremely high, in 
those few cases, it seems likely that while the model may be 
a good fit to the data it is not physically correct. 

In many cases, we find that the hot blackbody compo- 
nent contributes substantially to the overall flux. Substan- 
tial boundary layer emission is consistent with observed spins 
in neutron stars in LMXBs, which spin slower than the Ke- 
plerian angular speed in the inner accretion disk. If a neu- 
tron star had significantly higher angular speed, the bound- 
ary layer emission is expected to tend to zero as the angular 
speed approached the Kep lerian angular speed (see fig. 5 from 
iBhattacharvva et aT]|2000[ for example). 

6.5. Alternative emission line model 

We interpret the observed Fe K emission line as originat- 
ing in the inner accretion disk, with relativistic effects giv- 
ing rise to the asymmetric line profile (as is generally ac- 
cepted). There has, however, been a suggestion that these 
asymmetric Fe K lines can originate fro m an optically- 
thick high-velocity flow in these systems (Titarchu k et al.l 
2001 iLaming & Titarchufll2004t lLaurent & Titarchukl l2007l; 
Shaposhni kov et al.l I2009t iTitarchuk et al.l 120091) " In this 
model, a diverging wind with a wide opening angle is 
launched from the disk. The Fe K line is formed in a nar- 
row wind shell that is illuminated by X-rays from the accre- 
tion disk. The asymmetric profile arises when Fe K pho- 
tons scatter off electrons in the partly ionized wind, becom- 
ing red-shifted in the process. These photons undergo a large 
number of multiple scatterings to form the red wing of the 
line. Such a model appear s to fit some Fe K line profiles well 
dShaposhnikov et al.H2009l: iTitarchuk et al.ll2009l) . 



To reproduce the observed line profiles requires that there 
be a very large number of electron scatters i n the wind, and 
hence the wind had a high optical depth dTitarchuk et al.l 
120091 find T wind = 1.6 for neutron star LMXB Ser X-l and 
Twind = 4.9 for black hole X-ray binary GX 339—4). The mass 
outflow rate required for GX 339—4 is comparable to the Ed- 
dington mass inflow rate. In the low/hard state considered by 
Titarc huk et all d2009h . this means that the outflow rate is 30 
- 100 times higher than the inferred mass inflow rate. The 
outflow rate required to produce the line in Ser X-l is also 
comparable to the Eddington mass inflow rate. Not only are 
high outflow rates required, but the outflow must also not be 
coupled to the inflow rate given that we see similar line pro- 
files in the neutron star systems over more than two orders of 
magnitude in luminosity. 

Wide-angle, optically-thick outflows should have observ- 
able consequences in X-rays and in other bands. In X-rays, 
ionizing flux from the central engine should generate absorp- 
tion lines for any line of sight through the wind. Yet when 
absorption lines are observed in black holes, they are al- 
ways slow, with velocities less than approximately 1000 km/s 
(see, e.g., [M iller et al] l2006bl l2008al on GRO J 1655 -40 and 
lUedaet al.1 120091 on GRS 1915+105). In some cases, ab- 
sorption lines are seen in the absence of any relativistic line 
(e.g. H 1743-322 and 4U 1630-472. iMiller et al.ll2006cl 
iKubota et al.l 120071) . The absorption lines detected in most 
neutron star systems (the edge-on "dippin g" sources) are of- 
ten consistent with no velocity shift (e.g. iSidoli et al.l 120011: 
iParmar et alTM&fBomn et al.l2005tlbfaz Trigo et alj2006l) . 
Clearly, observed absorption spectra are in strong disagree- 
ment with the wind model. 

In optical and IR bands, optically-thick diverging outflows 
should also have consequences. Whether in a black hole, neu- 
tron star, or white dwarf system, the outflow should serve to 
obscure our view of the outer accretion disk and/or compan- 
ion star, at least during some intervals of the binary period. 
Along particular lines of sight, the outer disk should be ob- 
scured for nearly all of the binary period. Obscuration of 
this kind has not been reported, though it would surely be 
easily detected. Doppler tomography reveals no evidence of 
such outflows as accretion disk, stream/hot spot a nd compan- 
ion star emission are often clea rly revealed (e.g. Mar sh et al.l 
[19941 ISteeghs & Casaresll2002h . 

Finally, it is worth noting t hat relativistic lines are not de- 
tected in white dwarf systems. Tita rchuk et al.l d200 9) fit their 
model to the XMM-Newton CCD spectrum of GK Per, where 
they only observe one emission line. This, however, ignores 
the fact that high-resolution Chandra gratings observations 
clearly separated the emission into 3 separate Fe components , 
as is seen in many magnetic CVs dHellier & Mukai 20041) . 
The spectra of CVs is best described in terms of partially ob- 
scured emission from the boundary layer, which can be fit- 
ted wit h cooling-flow models that predict multiple iron charge 
states dHellier & Mukaill2004 . 

6.6. Equation of State 

To begin to put meaningful constraints on the neutron star 
equation of state, it is useful to combine both mass and ra- 
dius measurements. Of the sources in our sample, Cyg X-2 
is the only one with a well constrained m ass measurement, 
M NS = 1.5±O.3M dElebert et alJ 120091) . Combining this 
mass measurement with our stellar radius upper limit (8.1 
GM/c 2 ) can be seen in Fig. [10] This does not rule out any 
of the possible equations of state shown here, and tighter con- 
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TABLE 12 
Blackbody emitting radius 



Source 


Mission 


Obs. ID 


BB radius (km) 


BB radius (km) 








Phenomenological 


Reflection 


Serpens X-l 


Suzaku 


401048010 


5± 1 


5± 1 




XMM-Newton 


0084020401 


7±2 


7±2 




XMM-Newton 


0084020501 


6±1 


5±1 




XMM-Newton 


0084020601 


7±2 


6±1 


4U 1636-53 


XMM-Newton 


0303250201 


1.5±0.1 


1.0 ±0.2 




XMM-Newton 


0500350301 


4.0 ±0.1 


3.1 ±0.5 




XMM-Newton 


0500350401 


4.5 ±0.1 


3.8±0.3 


4U 1705-44 


Suzaku 


401046010 


3.5±0.4 


1 1+ 09 




Suzaku 


401046020 


5.2±0.2 


3.9 ±0.4 




Suzaku 


401046030 


3.9±0.2 


3.2±0.2 




XMM-Newton 


0402300201 


0.7 ±0.1 


0.7 ±0.1 


4U 1820-30 


Suzaku 


401047010 


5.4±0.3 


5.4±0.3 


GX 17+2 


Suzaku 


402050010 


3.1 ±0.1 


3.3 ±0.4 




Suzaku 


402050020 


4.2±0.2 


4.1 ±0.3 


GX 349+2 


Suzaku 


400003010 


12±3 


11 ±3 




Suzaku 


400003020 


10 ±3 


10 ±2 




XMM-Newton 


0506110101 


21 ±5 


19±5 


Cyg X-2 


Suzaku 


403063010 


3.1 ±0.6 


3.0 ±0.6 


SAX J 1808.4-3658 


Suzaku 


903003010 


9.0±0.8 


6.8±0.7 




XMM-Newton 


0560180601 


56±4 


540 ±74 



NOTE. — Assumed distances are given in TablefTI 



straints on mass and radius would be needed. 

Significantly improved sensitivity and spectral resolution 
from future X-ray missions, such as Astro-H and the Interna- 
tional X-ray Observatory (IXO) will certainly provide tighter 
constraints from modeling Fe K emission lines, and with the 
High Timing Resolution Spectrometer on IXO even the po- 
tential to study simultaneou s spectral and timi ng evolution on 
extremely short timescales (Barr et et al.ll2008h . 
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FIG. 10. — Constraints on the ultra-dense matter equation of state from the 
mass and radius upper limit of Cyg X-2 (light gray region). The dark gray 
region is not allowed due to causality. All mass-radius curves are as labeled in 
Lattimer & Prakash 12001) (see references therein for details of the equations 
of state ), except for NL4 which is from Akmal et al. 1 1998) and Z271 which 
is from Horowitz & Piekarewicz 12001). 
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